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PREFACE 


The Second Annual Report contained in this document covers the 
activities performed during the first two years of the NASA HOST 
Program, "Creep Fatigue Life Prediction for Engine Hot Section 
Materials (Isotropic)", under Contract NAS3-23288. The objective 
of this effort is to improve the high temperature crack initiation 
prediction technology for gas turbine hot section components. This 
program is being conducted under the direction of Dr. G.R. Halford 
who serves as the NASA Program Manager. The Program Manager and 
Principal Investigator at United Technologies Corporation (UTC) is 
Mr. Vito Moreno .J Mr D.M. Nissley significantly contributed to the 
analytical effort for data reduction and model evaluations. Dr. 
L.S. Lin performed all of the metal lographic examinations and 
interpretations of this information. The authors would like to 
express their gratitude to Messrs R. Masci and N. Glowac for their 
efforts in conducting the fatigue testing. 
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SECTION 1.0 


INTRODUCTION 

The overall operating cost of the modern gas turbine engine is greatly 
influenced by the durability of combustor and turbine structural components 
operating at high temperatures. Inadequate durability results in reduced 
engine efficiency and increased maintenance costs due to premature repair and 
replacement. To increase the durability of these components, more accurate 
structural analysis and life prediction methods must be developed for 
components operating at higher temperatures. However, improvements in the 
state-of-the art technology for elevated temperature durability prediction 
have been hampered by 1) the severe operating conditions of the engine, 2) the 
inability of analytical and life prediction tools used in the design of lower 
temperature components to predict complex material behavior and interaction of 
damage mechanisms of components at elevated temperatures , and 3) the high cost 
of engine development testing which prohibits the accumulation of adequate 
failure data and local operating conditions required for the systematic 
development and calibration of durability prediction models. 

Traditionally, component durability associated with fatigue cracking has been 
attributed to separate crack initiation and crack propagation processes. 
Current cyclic crack initiation prediction methodology is based on correlative 
fatigue models generally developed to address a specific component or loading 
configuration. These prediction models are directly calibrated to 
macroscopical ly observable and predictable quantities (strain range, mean 
stress, etc.) and do not address the specific mechanisms associated with the 
initiation process. This approach has been successful at lower temperatures 
but its applicability at higher temperatures, where interaction of deformation 
mechanisms and damage accumulation is strongly temperature dependent, has not 
been established. This is particularly true for the higher strength isotropic 
alloys used in gas turbine hot section structures. 


The NASA Hot Section Technology (HOST) program is aimed at developing improved 
life prediction technology for structures operating at elevated temperatures. 
As part of HOST, the present program will investigate fundamental approaches 
to high temperature crack initiation life prediction, identify modeling 
strategies and develop specific models for component relevant loading 
conditions. The program is a 5-year, 2-part effort (2-year base program plus a 
3 year optional program) that will consider two isotropic hot section 
materials and protective coating systems. Under the base program, various life 
prediction approaches for high temperature applications have been investigated 
and basic models for simple cycle, isothermal loading conditions selected and 
developed. Models that address thermomechanical cycling, multi -axial 
conditions, cumulative loading, environmental effects and cyclic mean stress 
will be developed under the optional program. These models will be verified on 
an alternate material. 
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SECTION 2.0 


SUMMARY 

This report summarizes the activities performed during the first two years of 
the NASA HOST Program, "Creep Fatigue Life Prediction for Engine Hot Section 
Materials (Isotropic)", being conducted by Pratt & Whitney Aircraft. The 
program is a 5-year, two part effort aimed at improving the high temperature 
crack initiation prediction technology for gas turbine hot section components. 
The two-year base program comprised the following tasks: 

Task I - Material /Coating/Component Selection and Acquisition 
Task II - Screen Candidate Life Prediction Approaches 
Task III - Evaluate Best Candidate Life Prediction Approach 
Task IV - Reporting 

Significant results of the program are listed below. A discussion of the 
technical progress of the program is provided in the following section. 

The major results of Task I are summarized below: 

1. Cast B1900 + Hf and wrought IN 718 were selected as the base and 
alternate materials, respectively. 

2. A single heat of B1900 + Hf was obtained and test specimens fabricated. 

3. The material was characterized with respect to grain size, 1 size, 
carbide distribution, and dislocation density. 

4. Monotonic tensile and creep testing demonstrated engineering properties 
within anticipated scatter for this material. 

5. Examination of the tensile tests has shown a transition from 
inhomogeneous "planar" slip within the grains at lower temperatures to 
more homogeneous matrix deformation at higher temperatures. 
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6. Examination of the creep tests has shown a transgranular failure mode 
at 760°C (1400°F) and an intergranular failure mode at 871 °C (1600°F) 
and 982 °C (1800°F). 

7. A study was conducted that investigated the effects of test specimen 
geometry and fabrication process on fatigue life. As a result, the 
axial strain controlled specimens were designed with a smooth (no 
extensometer ridges) gage section and fabricated using centerless 
grinding followed by light electropolishing. 

8. A fatigue test matrix was established to provide baseline data to 
define crack initiation life as a function of major variables and for 
life prediction model evaluation. Approximately 100 strain controlled 
fatigue tests have been completed. Major variables investigated were 
temperature, strain range, strain rate, mean strain, and strain hold 
times . 

9. Examination of specimens during testing indicated that measurable 
0.76mm(.030 in.) surface cracks appear early in the specimen cycle 
lives. This has been used as the definition of crack initiation for the 
initial model evaluation work. 

10. Observed crack initiation sites are all surface origins associated with 
either grain boundary carbides or local porosity. The initiation life 
is not significantly affected by the character of the site. 

11. Transgrannular cracking was observed at most initiation sites. 

The major results of Task II are summarized below: 

1. A total of seven different life prediction models were evaluated using 
the data generated in Task I. The specific models included: 

A. Coffin-Manson 

B. Damage Rate (Majumdar) 
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C. Hysteresis Energy (Ostergren) 

0. Crack Tip Opening Displacement 

E. Linear Elastic Fracture Mechanics 

F. Total Strain - Strainrange Partitioning 
u. Ductility Exhaustion 

Each model was evaluated with respect to its predictive capability and the 

types of data required for the determination of the model constants. 

2. Most of the models were able to accurately regress baseline data used 
for model constant determination. However, none were considered to 
demonstrate the degree of predictive capability together with simple 
data (test) requirements to serve as the basis for a practical and 
useable life prediction model for gas turbine hot section components. 

The major results of Task III are summarized below: 

1. A life prediction approach, designated Cyclic Damage Accumulation 
(CDA), has demonstrated good predictive capability with the minimum 
data requirements. As currently proposed all model constants are 
determinable from rapid cycle fatigue, monotonic tensile and creep 
tests. 

2. The model is fundamentally a Ductility Exhaustion approach in which the 
cyclic capability (ductility) of the material grains is exhausted by a 
damage parameter integrated through the cyclic evolutionary response. 

3. Tne transgrannular initiation, observed in the majority of tests, was 
considered to be indicative of a grain damaging process. Those cycles 
wnich generated intergrannular initiation were considered a grain 
boundary damaging process. 3ased on similarities of the deformation 
structures, the grain and grain boundary cyclic capabilities were 
associated with the amount of primary and secondary creep observed in 
monotonic creep loading tests. 
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4. At temperatures where time dependent deformation is normally not 
considered a significant process, the grain cyclic capability is 
assumed to be related to the amount of residual inelastic strain 
(% elongation) observed in a monotonic tensile tests. 

5. The cyclic damage parameter is based on determination of a reference 
damage rate. The reference condition for this work has been fully 
reversed, strain controlled fatigue tests conducted at a strain rate of 
0.00167 sec"^ . The reference damage rate is determined as the cyclic 
grain capability divided by the observed number of cycles to crack 
initiation. 

6. For cycles other than the reference condition, two damage ratios 
representing time independent and time dependent processes, are used to 
modify the damage reference rate. The damage ratios, determined from 
the stress response, are considered an advantage of this method since 
absolute damage increments are not calculated. 

7. Application of the model to the prediction of fatigue tests not used in 
the constant determination has shown the ability to predict strainrate, 
mean strain and hold time effects. 

8. Sensitivity analyses were conducted to evaluate the predictive 
capability for various assumptions regarding initiation crack 
definition and the accuracy of the local stresses. Some loss in 
predictive capability was demonstrated but the model was still able to 
adequately predict the trend in life associated with various cycle 
parameters. 
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SECTION 3.0 


TECHNICAL PROGRESS 

3.1 TASK I - MATERIAL/COATING/COMPONENT SELECTION AND ACQUISITION 

3.1.1 Hot Section Survey and Material Selection 

A survey was conducted of the isotropic materials and surface protection 
coating systems currently in use in the hot section of commercial gas turbine 
engines. The results of the survey, shown in Table I, indicate that the 
predominant hot section materials are high strength nickel base alloys in 
either cast or wrought form. Furthermore, the surface protection coating 
systems currently can be classified as either diffusion or overlay. A 
recommendation identifying the materials and coating systems for the base and 
option programs was submitted to the NASA Program Manager who approved the 
following selections: 

Base Material: Cast B1900 + Hf (PWA 1455) 

Alternate Material: Wrought IN 718 (AMS 5663) 

Coatings: Diffusion Aluminide (NiAl )/0verlay (MCrAlY) 

These selections provide an opportunity to investigate generic creep-fatigue 
crack initiation mechanisms and life prediction approaches for a wide range of 
relevant hot section materials and applications. The use of cast and wrought 
alloys allows the development and examination of life prediction models for 
two materials having the same matrix (Ni) but significantly different 
composition and microstructure (e.g., high vs low volume fraction V ). The 
selection of a diffusion and an overlay coating also provides a variation in 
composition and microstructure and will be useful in developing 
mechanistically based life prediction models. These material /coating systems 
provide the broadest range of generic information consistent with the goals of 
the HOST program. 
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Table I 


Hot Section Materials and Surface Protection Coating Systems 


Material Designation 

Form/Base 

Usage/Comments 

B1900 + Hf (PWA 1455) 

Cast - Ni 

o High Pressure Turbine (HPT) Blade 
and Vane Application 
o Potential Segmented Combustor Use 

IN 792 (PWA 1467) 

Cast + HIP - 

Ni o Si deplate and Vane Application 

IN 713C (PWA 655) 

Cast - Ni 

o Extensive Low Pressure Turbine 
(LPT) Vane Application 
o HPT Si deplate Application 

MAR-M-509 (PWA 647) 

Cast - Co 

o HPT 1st Vane Application 
o Strategic Material 

IN 718 (AMS 5662) 

Wrought - Ni 

o Extensive Projected Use in Hot 
Section Static and Rotating 
Structures 

Hasten oy X (PWA 1038) 

Wrought - Ni 

o Combustor Liner Sheet Alloy 
o Single Phase Alloy - Con- 
clusions may not be Relevant to 
Other Hot Section Superalloys 

Haynes 188 (PWA 1042) 

Wrought - Co 

o Combustor Liner Sheet Alloy 
o Strategic Material 

Coating Type 

Composition 

Comments 

Diffusion Aluminide 

NiAl 

o Extensive Blade and Vane 
Application 

o State of the Art for Oxidation 
Coatings 

Overlay 

MCrAlY 

o Different Composition and 
Microstructure 

o Compatible with Application of 
Thermal Barrier 
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3.1.2 Base Material Description 


The B1900 + Hf material selected for the program was part of a single heat, 
designated W-0098, obtained from Certified Alloy Products Inc., Long Beach, 
California. The chemical composition of this heat is compared to nominal 
specifications in Table II. A total of 2500 pounds of material was obtained 
for specimen fabrication. 


Table II 


Chemical Composition of B1900 + Hf (Heat W-0098) 

Element Nominal (%) Heat W-0098 


c 

0.11 

0.09 

Cr 

8.0 

7.72 

Co 

10.0 

9.91 

Mo 

6.0 

5.97 

A1 

6.0 

6.07 

Ta 

4.3 

4.21 

Ti 

1.0 

0.99 

B 

0.015 

0.016 

Zr 

0.08 

0.04 

Fe 

0.35* 

0.17 

w 

0.1* 

0.04 

Cb 

0.1* 

0.08 

Bi 

0.5 ppm 

0.1 

Pb 

10.0 ppm 

0.1 

Hf 

1.5 

1.19 

Ni 

Remainder 

Remainder 


♦Maximum 

Two casting geometries were selected for specimen fabrication. The "standard" 
bar shown in Figure 1A was initially selected because it could accommodate all 
anticipated test specimen configurations. The constant section bar (Figure IB) 
was later cast to investigate the effects of porosity on material properties 
(see Section 3.1.4). For both bar geometries the casting parameters, melt and 
mold temperatures, were established to produce a uniform small grain size in 
each bar. This would result in a large number of grains in a specimen 
cross-section and insure an "isotropic" material response. 
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30.5 
<1.20)DIA 



(A) "Standard" Bar Configuration 



(B) Constant Section Bar 


Figure 1 Casting Geometries Used for Test Specimens 

All bars were fully heat treated prior to machining. The heat treatment cycle 
included: 

Solution - 1079 + 14°C (1975 + 25°F) for 4 hours; air cool 

Precipitation - 899 +_ 14°C (1650 + 25°F) for 10 hours; air cool 

For future reference, all cast bars were given an identifying designation. The 
"standard" bars were cast four bars at a time. The molds were numbered 
sequentially from number one, with individual bars designated A, B, C or D. 

The cylindrical bars were numbered sequentially starting with 647. 
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The structure of the material was documented in both the as-cast and fully 
heat treated condition using optical, SEM and TEM techniques. The following 

observations were made: 


1 . 


Use of a micrograph (300X) showed the grain size of both the as-cast 
and fully heat treated material to be between ASTM No. 1 and 2. This 
represents an average grain diameter of 0.018 cm (0.007 in) to 0.025 cm 
(0.010 in) which should produce isotropic stress and strain conditions 
in the test specimens. A comparison of the micrographs and the ASTM 
standard is shown in Figure 2. 


2. The replica technique was used to study the gamma prime, Y 1 , size and 
distribution within the grains. In the as-cast material, the Y size 
was 0.6 nm while the fully heat treated material showed an increase in 

size to 0.9 pun (see Figure 3). 


3 MC carbides were observed in the grain boundaries of both the as-cast 
and fully heat treated material (see Figure 4). Microprobe analysis 
indicated that the carbides are enriched in Ta, Ti, Mo, and Hf. 


4. For future assessment of deformation accumulated in specimens after 

testing, initial dislocation density measurements were made. 

Representative dislocation networks for the as-cast and fully heat 
treated materials are shown in Figure 5. Calculated densities^ about 
9.3 x 10 4 * * * * 9 /cm 2 for the as-cast material and about 5.7 x 10 /cm 
for the fully heat treated material. 


3.1.3 Tensile and Creep Characterization 

Honotonic tensile and creep testing of B1900 + Hf specimens was conducted to 

document typical engineering quantities and define deformation and failure 
mechanisms. Specimen geometries used for this testing are shown In Figure 6. 
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As-Cast y'= 0.6 ^m 



Figure 3 Gamma Prime (V ) Shows Increase with Heat Treatment 
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Fully Heat Treated 

Figure 4 Typical Distribution of Carbides Along Grain Boundaries 
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Fu77y Heat Treated P= 5.7 x 70^/cm^ • 


Figure 5 Dislocation Networks for As-Cast and Fully Heat Treated Materials 
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Figure 6 Test Specimens for Monotonic Tensile and Creep Testing 
Tensile Tests 


A total of 21 monotonic tensile tests were conducted. A summary of all test 
conditions and observed properties is presented in Table III. A majority of 
the tests were run at the ASTM standard rate of 0.005 min. -1 . Several tests 
were also conducted at 0.0005 min"^ to investigate rate effects on the 
tensile properties. As indicated, a number of tests terminated with fracture 
occurring outside of the gage section of the specimen. Despite these failures, 
most of the properties measured early in these tests (e.g., modulus, 0.2% 
yield) are consistent with results obtained from tests in which failure 
occurred within the specimen gage section. An examination of these specimens 
suggested that the tendency to fail outside the gage section is associated 
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Table III 


Summary of Tensile Testing 


Temp. 

:(°f) 

Spec. 

ID 

min 1 

ExlO" 3 

MPa(KSI) 

.2% Yield 
MPa(KSI) 

Ult. 

MPa(KSI 

Elong 

% 

RA Comments 

% 

r 

18A 

0.005 

187.5(27.2) 

714(103.5) 

- 

4. 1 3 

5.9 3 

50(500) 

18B 

0.005 

169.6(24.6) 

702(104.8) 

888(128.8) 

8.3 

10.7 


38(1000) 2A 0.005 149.6(21.7) 727 ( 105j 

0.005 143.4(20.8) 701(11 


49(1200) 19A 

60(1400) 2C 0.005 

2D 0.005 

38B 0.005 

650 2 0.005 

38A 0.0005 

647 0.0005 


(A 


150.3(21.8) 709(] 

122(17.7) 721 

144.1(20.9) 739i 

157.8(22.9) 70 
122.7(17.8) 77 
146.8(21.3) 7Q 


4.7 7.2 

^6(115.5) - 

1(138.5) 8.0 9.1 

947(137.4) 7.8 7.6 


71(1600) 19B 0.005 146.9(21.3) fl 

39B 0.005 

651 0.005 

39A 0.0005 

648 0.0005 

82(1800 ) 20A 0.005 

41B 0.005 

652 0.005 

41A 0.0005 

649 0.0005 

1093(2000) 20B 0.005 


147.5(21.4) 

131(19.0) 

156.5(22.7); 

W 

KTC m 

792(114.8) 

5.7 

6.8 

146. 8(21. 3M 

mm 

780(113.1) 

5.7 

5.3 

111(16. iM 
135. 8 (IK 
119. 3(A . 

’ jH®> 

iWm.6) 

. 3 ) 
. 7 

470(68.1) 

492(71.4) 

6.5 

7.7 

5.0 

8.8 

144 .if J 

v * * * 
4 1.1) 

471(68.3) 

6.5 

7.2 

101 .K* 

91(24.3) 

- 

- 

- 


| fractured outside gage sectj 
- specimens 648 - 652 machinar 
j within representative scatj 


0 0 

Co 


inch diameter bars 


stopped at 
1.5% for ex 


0.G J 

O.G. 

O.G. 


O.G. 


stopped at 
for exam. 
O.G. 

stopped at , 
for exam. 


stopped at 1 
for exam. 


stopped at 1 
for exam. 


stopped at 1 
for exam. 


with a higher degree of porosity at the thicker ends of each casting as 
illustrated in Figure 7. To further investigate these failures, a series of 
specimens was fabricated from 15.9mm(5/8" ) diameter cylindrical bars. While 
cast with the same parameters as the previous castings, the cylindrical shape 
was thought to produce a more uniform distribution of porosity. Six tensile 
tests were run with these specimens (designated 647 to 652). As shown in the 
table, the tensile properties are consistent with previous tests and all 
failures occurred in the specimen gage section. 

Comparisons of the 0.2% yield strength and measured static modulus of 
elasticity vs the anticipated representative scatter in B1900 + Hf are 
presented in Figures 8 and 9. Representative stress - strain responses at 
various test temperatures are shown in Figure 10. 
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Figure 7 Observed Area of Higher Porosity in As-Cast Bar 



" — 0.2% YIELD STRESS 
O "STANDARD CASTING" 

□ 1 5 9 MM IB/0 INI DIA BAR CASTING 


097 (1301 
028 (120) 
759 (110) 

690 (100) 

E75 

* 621 (90) 

< 


» B Q \ 

19A O \ 





Figure 10 Monotonic Tensile Response - e» .005 min 


Six specimen tests were terminated at approximately 1 - 1.5% strain for 
examination to characterize the change in dislocation structure with 
increasing temperature. As expected, the lower temperature < 871°C ( < 1 600°F) 
deformation is characterized by Inhomogeneous slip as Illustrated by the 760 # C 
(1400°F) deformation shown in Figure 11. At higher temperatures , as transition 
between the Inhomogeneous "planar* slip to more homogeneous matrix deformation 
in the grains is observed. The deformation at 1093°C (2000°F) is shown in 
Figure 12. A composite for all the temperatures examined at the same strain is 
presented in Figure 13. Note that microtwins were observed in the latter 
stages of deformation (=“4.5%) at 760°C (1400°F). 


20 



ORIGINAL P^V“j:; r 
OF. POOR QUALi i 4 



Li_ 

o 

O 

O 


rO 


C 

o 


$- 

o 


<L> 

Q 

S- 

C 


a> 

s- 

3 

CD 

•r* 

U- 



21 


Figure 12 Homogeneous Matrix Deformation at 2000°F 
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Figure 13 Composite TEM Photos of Te 
Various Temperatures 


A. 538°C (1000°F), 1.5% STRAIN 

B. 760°C (1 400°F), 1.3% STRAIN 
760° C (1 400°F), 4.5% STRAIN 

C. 871 °C (1 600°F), 1 .3% STRAIN 

D. 982°C (1 800°F), 1.2% STRAIN 

E. 1093°C (2000°F), 1.0% STRAIN 


TENSILE 


sile Tested B- 7 900 Hf Specimens at 



Creep Tests 


A summary of the test conditions and observed properties for all monotonic 
creep tests conducted is presented in Table IV. As indicated, the tests were 
conducted at four temperatures [649 V C (1200T), 760“C (1400D, 871 C 
(16UOD, 982 W C (1 BOO'D] and at stress levels representing 50 to 122% of the 
bottom of the yield stress scatter band. As shown in Figure 14 the normalized 
stress {< 7 /a 2 % Ys ) reduces the apparent temperature effect for specimens 
having the same deformation mode, i.e., specimens tested at 871 C (1600 F) and 
980 V C (1800D have similar rupture lives for the same normalized stress and 
intergrannular failure, while specimens tested at 760 C (1400 F) and 649 C 
(1200D show a significantly longer life for the same, or higher, normalized 
stress level. 

Table IV 

Summary of Creep Test Results 

Secondary 


Spec- 

Temp. 

Stress 

% Min. 

Frac- 

Life 

Creep Rate Elong. 

RA 

imen 

-cm 

MPa(KSI) 

Yield 

ture 

(Hrs. ) 

(Min-' ) 

% 

% 

19C 

982(1800) 

234(34) 

75 

I 

20.6 

2.5 E-05 

6.0 

7.8 

190 

982(1800) 

283(41 ) 

90 

I 

4.1 

0.7 E-04 

3.0 

2.4 

39C 

982(1800) 

283(41) 

90 

n.a. 

2.9 

2.0 E-04 

5.5 

5.5 

41 C 

871(1600) 

283(41 ) 

50 

n.a. 

441 

6.0 E-07 

2.6 

n.a. 

18C 

871(1600) 

427(62) 

75 

I 

18.2 

2.5 E-05 

3.2 

4.9 

40B 

871(1600) 

427(62) 

75 

n.a. 

20.3 

1.5 E-05 

2.26 

2.7 

18D 

871(1600) 

517(75) 

90 

I 

2.8 

1.0 E-04 

3.2 

5.4 

40A 

871(1600) 

517(75) 

90 

n.a. 

2.5 

1.25 E-04 

2.3 

3.5 

598 

871(1600) 

586(85) 

104 

n.a. 

0.68 

3.85 E-04 

2.87 

5.04 

47D 

871(1600) 

586(85) 

104 

n.a. 

0.44 

7.80 E-04 

2.79 

n.a. 

1 

871(1600) 

586(65) 

104 

n.a. 

0.24 

1.40 E-03 

2.72 

n.a. 

7B 

871(1600) 

690(100) 

122 

n.a. 

0.082 

6.40 E-03 

3.00 

2.40 

7C 

871(1600) 

690(100) 

122 

n.a. 

0.047 

1.54 E-02 

3.22 

n.a. 

408 

871(1600) 

690(100) 

122 

n.a. 

0.035 

1.25 E-02 

1.97 

n.a. 

697 

760(1400) 

600(87) 

90 

n.a. 

38.5 

2.64 E-06 

1.47 

2.41 

2UC 

760(1400) 

600(87) 

90 

T 

134.8 

n.a. 


3.1 

200 

760(1400) 

669(97) 

100 

T 

30.2 

7.9 E-06 

2.9 

9.4 

390 

760(1400) 

669(97) 

100 

n.a. 

49.8 

7.5 E-05 

3.78 

5.9 

o02 

649(1200) 

690(100) 

102 

n.a. 

1958.7 

1.18 E-07 

1.10 

1.60 


I = Intergranular 
n.a. = Not Available 
T = Transgranular 


Examination has shown that at the higher temperatures [871 °C (1600°F) and 
982 °C (1800 # F)] the specimens failed In an Intergranular cracking mode, while 
the specimens tested at 760°C (1400°F) showed a greater tendency for 
transgranular cracking (see Figures 15 through 18). Substantial differences In 
the grain dislocation structure were also observed. In specimens tested at 
932°C (1800°F) and 871°C (1600°F) (Figures 19 and 20), dislocation networks 
surrounding the gamma prime (y‘) particles are the dominant features. This Is 
an indication of a loss of coherency of the r' and the matrix. At 760°F 
(1400°F), no dislocation network is observed; Instead short segments of 
dislocation and microtwins (similar to the tensile tests) are inhomogeneously 
dispersed throughout the grains. The presence of the microtwins at this 
temperature may account for the tendency of cracks to progress in a 
transgranular manner. 

A comparison of the percent elongation observed in several of the specimens 
during the monotonic tensile and creep testing Is presented in Figure 21. As 
indicated, the creep elongation is generally smaller than the tensile 
elongation. 



Figure 14 Variation of Creep Rupture Life with Stress 
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CREEP 

FRACTURE MODE 


982°C/283 MPA 
(1 800°F/41 KSI) 



871°C/517 MPA 
(1 600°F/75 KSI) 


INTERGRANULAR 



TRANSGRANULAR 


760°C/600 MPA 
(1400°F/87 KSI) 


Figure 15 Creep Fracture Modes 
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section along line AA' showing intergrannular fracture by SEM. 
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surface (b) Longitudinal section along line AA' showing 
intergrannular fracture by optical microscopy (c) Longitudinal 
section along line AA' showing no indication of intergrannular 
fracture by SEM. 
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Figure 21 Elongation Observed in Tensile and Creep Tests 


3.1.4 Fatigue Specimen Design, Fabrication and Test Facility 

The three axial, strain control specimen designs used in the initial fatigue 
testing are shown in Figure 22. Type A has integral, large extensometer 
ridges; type B has integral "mini" ridges; and type C is a smooth cylindrical 
gage section design. An initial group of 18 specimens was fabricated to 
investigate machining and geometry effects on fatigue life. Six type A 
specimens were electrochemically ground and six each of types B and C were 
centerless ground from "standard" B1900 + Hf castings (see Figure 1 A) . All 
specimens were then electrochemically polished. Fully reversed, strain 
controlled fatigue tests were conducted at: 

871 °C (1 600°F ) , Ae= + 0.25%, e = 1.67 E-03 S" 1 

538°C (1000°F) , Ae= + 0.25%, e= 1.67 E-03 S' 1 
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Both total specimen cyclic life and failure location were recorded. Two 
specimens were subsequently eliminated from the stucly because of testing 
problems. A Ueibull analysis of the data, presented in Figure 23, shows that 
the 871 °C (1600°F) data for each specimen type is uniformly distributed, 
indicating little effect of geometry or machining procedure at this 
temperature. Primary failure locations occurred throughout the gage section 
for all specimens tested. However, at 538°C (1000 # F), the type C data are 
grouped at the higher end of the distribution. Furthermore, the type A and B 
failure locations were predominantly associated with the extensometer ridge 
fillets, whereas the type C failures generally occurred throughout the gage 
section. These results suggest a fillet or "notch" effect in the fatigue life 
at the lower temperature. The distribution of the type A and B tests also 
suggests that the machining procedure (ECG vs CG) has little effect on fatigue 
life at these conditions. 

S38°C 



Figure 23 Weibull Analysis of Initial Fatigue Data to Investigate Geometry 
and Machining Effects 
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Further investigation using X-ray diffraction to measure residual surface 
stresses produced by both processes resulted in similar average stress levels, 
but significant amounts of scatter were observed in the measurements. The 
X-ray diffraction (XRD) technique is successfully used on fine grain (wrought) 
materials; however, application to a cast material, which has a relatively 
larger grain size, results in fewer grains being sampled by the incident X-ray 
beam, and thus, much more scatter in the stress measurements. Data taken on 
several B1900 + Hf specimens machined using centerless grinding and 
electrochemical grinding is shown in Figure 24. For comparison, measurements 
taken on fine grain (wrought) IN 718 material show considerably less scatter, 
with centerless grinding producing a slightly larger compressive stress. 

A second group of specimens was fabricated with the type B geometry using 
centerless grinding and a light [< .025mm(<0.001 in.) removal on a diameter] 
electropolish. The cyclic lives of these specimens, tested at the same 
conditions, showed similar trends to the previous data. Subsequent specimens 
were fabricated to the type C geometry in order to eliminate the notch effect 
at lower temperatures using centerless grinding and a light electropolish. 

CG - CENTERLESS GRINDING 

ECG - ELECTROCHEMICAL GRINDING 

CAST 61 900 + Hf WROUGHT IN7 18 


RESIDUAL CG ECG H4718 IN718 



Figure 24 Residual Stresses Measured by X-Ray Diffraction 
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The test facility used for the isothermal fatigue tests included a servo- 
controlled, closed loop hydraulic testing machine with MTS controllers, a 7.5 
kw - lOkhz Tocco induction heater, and an Ircon infrared radiation pyrometer 
for temperature measurement. Induction heating was selected to facilitate 
setup and inspection of the many fatigue tests. The Ircon pyrometer is used to 
measure temperatures above 760°C ( 1 400 °F ) ; below this temperature, a 
thermocouple monitored specimen temperature. Axial strain measurement was 
accomplished with an MTS extensometer. The quartz rods which define a one-inch 
gage section are spring loaded against the specimen and did not show any signs 
of slippage during testing. A typical test setup is illustrated in Figure 25. 



Figure 25 Typical Test Setup for Isothermal Fatigue Testing 
3.1.5 Baseline Fatigue Testing 

Isothermal, strain controlled fatigue tests were conducted to define the crack 
initiation life of B1900 + Hf material and to provide a baseline for life 
prediction model evaluation. Testing was limited to key variables relevant to 
the general creep fatigue life prediction problem. A general matrix of the 
tests is presented in Figure 26 and includes strain range (Ac), strain rate, 

ORIGINAL PAfc;^ ^ 

OR POOR QUALH i 


strain ratio (Re), dwell periods and temperature. All of the tests conducted 
for the base program are summarized in Appendix I. Those tests conducted at a 
strain rate of 1.67 E-03 sec ' represent the upper limit of rates tested and 
are considered to be fast enough to preclude a significant creep fatigue 
effect. Representative hysteresis loops from 871 °C (1600°F) tests conducted at 
strain ranges of 0.5% and 0.8% indicate a small amount of inelastic strain 
relative to the total strain range (Figures 27 and 28). A review of the cyclic 
response histories of the Rc = -l tests indicates that at 871 °C (1600°F) a 
small amount of cyclic softening occurs, while at 538°C (1000°F) cyclic 
hardening is observed at the larger strain ranges. The representative 
inelastic strain ratio of four specimens and the range of values for all 
specimens tested at the same conditions are presented in Figures 29 and 30. 


TEST 

CONDITION 

1 

2 

3 

STRAIN RANGE 
Az 1 Ac 2 4*3 
X 

X 

X 

STRAIN RATIO 

— oo _ 1 0 

X 

X 
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STRAIN RATE 

«1 h ‘3 

X 

X 

X 

DWELL PERIODS 

‘MAX ‘MIN ‘INT 

1l t 2 lj 1 2 t, t 2 

TEMPERATURES 
T 1 T 3 T 4 T 5 

X 

X 

X 
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X 

X 

X 


X 

5 

X 

X 
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X 
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Figure 26 Test Matrix for Baseline Fatigue Tests 
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MP* 414 CYCLE NO. 



Figure 27 


Cyclic Response 


of 871°C (1600°F) Fatigue Test with Ae* + .25* 



Figure 28 Cyclic Response of 871°C (1600°F) Fatigue Test with A« » + .4* 





± 0 . 25 % 



Figure 29 Inelastic Strain Ratio for 871°C (1600°F) Fatigue Tests 
(e = 1.67 E-03 S" 1 ) 



a 



0 .2 .4 .6 .8 1.0 

LIFE FRACTION (n/N) 


Figure 30 Inelastic Strain Ratio for 538°C (1000°F) Fatigue Tests 
( c* 1.67 E-03 S" 1 ) 


For each specimen test, the crack Initiation life was defined by three 
methods; (1) total separation of the specimen, (2) cycles to 5% and 10% stress 
amplitude drop from the steady state values, and (3) first Indication of 
cracking. 
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The first indication of cracking was established based on inspection of a 
number of specimens during testing. At each inspection, an acetate replica of 
the gage section was taken and examined to determine the number and size of 
surface cracks. In those tests where multiple cracks were observed, the crack 
that ultimately led to specimen failure was tracked throughout the life. At 
least one specimen was inspected for each condition shown in Appendix I. A 
summary of the fully reversed data based on the longest observed crack at each 
inspection is presented in Figure 31. General observations include the fact 
that measurable cracks 0.13 - 0.26mm (.005 - .010-in) appear earlier (based on 
percentage of separation life) and have faster growth rates at the higher 
temperatures. Specific details associated with each test type are discussed in 
Section 3.3. 

During the initial life model evaluations, crack initiation was defined as the 
generation of a 0.76mm (0.030-in) surface crack. Inspection of several fatigue 
fracture surfaces has shown that the assumption of a semicircular crack to 
estimate crack depth is reasonable for surface crack lengths to « 2.5mm 
(»0.100 in.). Thus, the 0.76mm (0. 030-in) 'surface length represents a 0.38mm 
(0.015 in) depth crack as compared to a material grain size of 0.18 - 0.25mm 
(0.007 - 0.010 in) (see Section 3.2.1). 

The original eighteen fatigue specimens, used to Investigate geometry and 
machining effects on fatigue life, were also examined to determine the 
location of the primary failure initiation sites and the mode of cracking. All 
specimens were found to have surface initiated fatigue cracks associated with 
either carbides and/or porosity. Two representative specimens and failure 
sites, as identified by SEM and microprobe analysis, are presented in Figures 
32 and 33. Failure sites determined to have a high tantalum (Ta)content or 
other constituents found in the grain boundary carbides were identified as 
carbide related failures. Failure sites having a chemical composition similar 
to the matrix, and where a pore could be observed, were identified as porosity 
related failures. Analysis of the specimen life data (Figure 34) suggested 
that the nature of the Initiation site (carbide vs porosity) is not 
statistically significant in determining the fatigue life. Based on this 
result, carbide or porosity initiation was not considered a primary variable 
in the analysis of subsequent fatigue tests. 


40 



SURFACE CRACK LENGTH, IN SURFACE CRACK LENGTH, IN 




ORIGINAL PAGS & 
OF POOR QUALITY 










CO 


CD 






h- 


C 





• 



•r— 


O 



!S 

0- 


<d 


CO 



O 

o 


CD 


s- 



s- 




+1 

o 


u 

s- 

c 



** 



co 

to 

o 


CD 

LO 

o 


s- 


*f— 


4-> 

C\J 

3C 


u 

>> 

-M 


4- 

• 

+-> 



-Q 

3 


C 

o 



>> 


X) 




CD 


S- 

a 

♦r- 

. . 

Q 

+ 1 

c 


to 

* 

s- 

CQ 

LD 


•r— 


E 

o 

-M 



II 



•r— 

JZ 

00 


C 

4-> 

o 

CO 

S- 

00 

•f— 

O 

CD 

w 

-C 


Q_ 


"O 

+-> 

E 


to 

“O 


00 



•r- 



c 

+-> 

*r— 

CD 

CD 

CJ 

->■ N 

00 

CO 

CO 


-C 

C 

CD 

U- 





+-> 

•r— 

Q. 

o 

• CL 

c 

CD 

C 


"O 

t/> 

O oo <0 


O 


CD 

c 


O CD i- 

00 

CO 


C 

o 

<4- 

LD ■— O 

CD 

0- 

CD 

•r* 

Q. 

O 

»— o O 

4-> 

s- 

+-> 


00 


'* — - ' 

>> +-» 

•r— 

3 

*r— 

o 

CD 

CD 

o o 

00 

00 

00 

J= 

U 

O 


CO 




00 

J- 

CO 

o lo i- 

C 

c 

C 


o 

<4- 

»— LO *4- 

o 

o 

o 

CL 

u 

t- 

cd 

•1“ 


•f— 

CO 


3 

00 C\J ■ — 

+-> 

CD 

-M 

E 

CO 

t/> 


<o 

CO 

T3 

cO 


CD 


S- O 

*i— 

•r— 

•r— 

>> 

s- 

Of 

CO O *r- 


-Q 

-M 

<o 

cO 


<4- CL 

•r— 

t- 

*r— 

u 


3 

■o 

>0 

C 

CO 

C 

1 

c 

+-> 

CD E 1— 

•r— 

C_> 

•r» 

X 

•f— 

O 

-M CL 






fO 

00 o — - 


<—>. 




S- 

<D O CO 


JD 


CJ 


U- 

h- «— 


** — 


> — 



CSJ 

CO 

< 1 ) 

S- 

3 

o 

•r* 



42 




OREGKSSAL PAL.L £$ 
OE POOR QUALITY 



CD 





C 




c 





o 




•1— 





•r— 




a> 





-4-> 




CO 





ro 


-C 



•* 




•r— 


4-> 


S- 



u 




•r— 


a> 

LD 


<T3 


•f— 


2 


+-» 

C\J 


s- 


C 




4 - 

• 


o 


•r- 


co 


*£ 

o 






c 



+ 1 


CO 




o 


Q 


c 


U 


r — 




•f— 


rO 


ro 



II 


2: 






c 

4-> 


o 


u 


CO 


cd 



-C 




c 


H 

<] 


CO 


-M 


o 


•r— 





ro 


•r- 


u 

, — V 

* 

to 

<D 



4-> 


cd 

Ll 

CO 

c 

O 

CD 


ro 


CL 

0 

CD 

CL 

ro 

CJ 


*r— 


LTi 

o 

r— 

<T5 

4- 

rO 





o 

O 

t. 

S- 

4— 


-M 

to 

4- 

o 

>> 

CO 

13 



CO 


O 


CJ> 

o 

to 

Z3 



u 


v — " 


4-> 


CO 


CD 

rO 

CD 


o 

o 

c 




i- 

O 

C_J> 

CD 

ro 

o 

c 


cn 

o 

<o 

o 

LO 

5- 


o 


*r— 


4— 

00 

* 

4- 

"O 



+-* 

>> 

S~ 

CO 

C\J 


CD 

>> 


rO 


13 

LD 

r— 

r— 

+-> 

+-> 


4- 

ro 

00 



ro 

rO 

•r- 



"a 


+J 


O 

*r- 

to 


4-> 

c 

CD 

(O 

o 

•r— 

-M 

o 

CD 

C 

o 



4— 

CL 

*r— 

s- 

-M 

•i — 

o 


"O 


>> 

C 

o 

•r— 

rO 

CD 

+-> 

CD 

E 

1— 

•r— 

a. 

CO 

Ll 

to 

o 

-M 

CL 







ro 

CO 

o 

„ — , 


^ — „ 


„ — ^ 


u 

CL) 

o 

»T3 


X> 


O 


Ll. 

1— 

f-~ 

v — 


’ — " 




CO 









CO 









CD 









t- 


















co 









•r— 









u . 











43 





Figure 34 Distribution of Initiation Sites for Initial Fatigue Tests 


In order to clearly identify the process of crack growth from the Initiation 
site to the 0.38mm(0.015 in.) crack depth defined as initiation, several 
specimens were ground down longitudinally until the Initiation site was exposed 
in the plane of polishing. Figure 35 shows the longitudinal view of a carbide 
initiation site in specimen 3D tested at 538°C (1000°F). The propagation path 
appears to be transgranular for approximately 4.8mm(0.190 in) (~20 grains) 
before switching to an Intergranular mode. In Figure 36, examination of 
specimen 6B tested at 871 °C (1600°F) identified the initiation site as a 
0.005mm ( 0.0002 in) surface pore. The crack grows approximately 0.76mm (0.030 
in) (3 grains) transgranularly before switching to an intergranular mode. 
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Examination of additional specimens has shown an increasing tendency toward 
intergranular cracking with higher temperatures and slower strain rates. 

Figure 37 shows the representative crack growth path in specimens tested at 
538°C (1000°F) and 871 °C (1600°F). The approximate crack depth for 10% and 50% 
of the specimen test life is superimposed. Clearly, a large portion of the 
cyclic life is spent on transgranular crack growth at these conditions. 

Transmission Electron Microscopy (TEM) was used to characterize the 
dislocation structure developed by the cyclic testing. The typical dislocation 
structure of a specimen tested at 538°C (1000°F) (4D) is shown in Figure 38. 

It appears that the large Y' particles impede the motion of dislocations which 
results in the dislocations tightly encasing the individual particles of Y'. 
Two active octahedral slip systems were identified: (111) <110> , location A, 
and (111) <101> , location B. The typical dislocation structure of a specimen 
tested at 871 °C (1600°F) (6D) is shown in Figure 39. In this case the 
dislocations do not form the tight encasement of Y' particles; rather, they 
gather in the matrix regions along the direction of the Burgers vector. 

Since cross slip is easier at 871 °C (1600°F) than at 538°C (1000°F), as the 
dislocations encounter resistance from the Y 1 particles they apparently 
overcome the resistance by cross slipping onto a different slip system. There 
appears to be a general directional alignment of the dislocation activity in 
Specimen 6D (see Figure 39). Three active octahedral slip systems were 
identified in this specimen: (111 ) < 1 1 0 > , (111 ) < Oil > and (111 ) <101 > , 
typical of dislocations at locations A, B and C. 

The thin foil specimens of 4D and 6D were taken ~7.1mm(~l/4 in) below the 
fracture. In order to establish that the dislocation structures observed were 
representative of deformation at the fracture surface, a thin foil was made on 
a second specimen tested at 538°C (1000°F), specimen 3D, Immediately below the 
fracture. The results (Figure 40) generally agree with the observations made 
on Specimen 4D which suggests that the deformation process is at least 
initially "homogenous" throughout the gage section. Specific observations used 
in the development of the proposed life prediction model will be discussed in 
Section 3.3. 
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crack initiation site is shown by arrow, (b) Crack propagates 
transgrannularly and becomes intergrannular at location shown by 
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Figure 37 Primary Crack Depth at 10% and 50% of Specimen Cyclic Life 
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Figure 38 Dislocation Structures After 
Being LCF Tested at 538°C( 1000°F) , lOcpm 
A«t = + 0.25% for 12,500 Cycles, Specimen 4D 
(~ 1 / 4 in. away from fracture surface). 

Zone axis (211) i A: (111) <lTo> 

\ B: (111) <101 > 


i9 



ORIGINAL PAGE 'IS' 

OF POOR QUALITY 


4-5 


n3 



LO 

■o 

VO 

<D 

Ol 

4-5 

•k 

00 

CNJ 

0) 


1— 

fc. 


o 

u. 

*4- 

CJ> 


— J 

K 

CD CNJ i 

c 

* 

•r— 

O J 

<U 


CO 

+ 1 i 


II 

<L> 

4-> 

4-> ! 

0- 

< 

<* 

<1) 

W* Vi 

t- 

E 

3 

1 

4-» 

U 

8 ' 

3 

r— 



4-5 

* 

i/1 



U. 

c 

0 

o 

O 

•i— 

° t 

4-> 

VO 

<0 

r— • 

o 


o 

, 


O r 


r— j 

•r- 

r-*. 

Q 

00 t 


Cycles, Specimen 6D (—1/4 in away from 
fracture jurfage). 

A: (111) <110>j 


tW) 


0R5C.-- ' 

OF. POOR QOfwO 



S8 

<8 "oJ 

. — jo 


Ll. C 
o •*— 



^ o 


E 


CO E 
CO CM 
LO LO 


<T3 > — I 


"O Q 
<D CO 


00 c 



U_ O 
O <1> 
_J CL 

oo 

CD 


•r- 00 
<D CD 
CO •— 

s- >> 

<D O 
+J 

4- O 
< CO 
CO 
to * 
CD CO 


<J O 


n lo 
CSJ 
c • 
o o 


fO 

U II 

o 

r- -M 

00 yy 

Q <3 


o 


<D 


cn 


51 


fracture surface]. 


3.2 TASK II SCREENING OF CANDIDATE LIFE PREDICTION APPROACHES 

3.2.1 Model Selection and Screening Criteria 

One of the fundamental objectives of this research effort is the 
identification of a fatigue life methodology or approach suitable for the wide 
range of conditions experienced by hot section structural components. For this 
effort, it is necessary to distinguish between an “approach" and a "model". A 
life prediction approach is considered as a conceptual basis for the 
prediction of cyclic life. Within current life prediction technology, 
approaches span the range of purely phenomenological to mechanistic. A life 
prediction model is the application of an approach in a mathematical form (a 
function or functional) for the prediction of life. On a relative basis, 
existing models can be characterized as shown in Figure 41, where the extreme 
left represents a purely phenomenological approach with increasing degree of 
mechanistic basis associated with the position to the right. As suggested in 
the figure, a model such as the Frequency Modified Life (Refs. 1,2), which 
relates the fatigue to inelastic strain range and frequency, is empirically 
derived using the observable macroscopic variables. All of the local 
macroscopic effects which determine life (crack initiation) are embodied in 
these parameters which are considered constant through the cyclic history. 



Figure 41 Characterization of Life Prediction Models 
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The Strainrange Partitioning (SRP) (Refs. 3,4,5) model also relates cyclic 
life macroscopic (observable) parameters, but has, as its foundation, the 
notion of local damage associated with reversible rate dependent effects. For 
this reason, it is considered to be more mechanistic than Frequency Modified 
Life. Chaboche's Damage Rate (Refs. 6,7,8) model represents a different 
approach from the previous two models. Local damage is associated with creep 
rupture and rapid cycle fatigue which accumulate non-1 inearly. Finally, the 
model of Hajumdar (Refs. 9,10,11) relates the local damage to microscopically 
observable details (cavity growth, grain boundary cracking, ...). Here, the 
life relationships are represented as functions of the local damage process. 
This model is considered more mechanistically based than the previous three 
model s. 

From a practical standpoint all of the models use macroscopically observable 
quantities determined from tests. The primary difference in the approaches is 
the degree of mechanistic reasoning or justification. None of the models 
considered use explicitly measured micro-mechanistic properties. 

During the initial screening, models representing the entire range of 
approaches were considered for selection of the best approach for prediction 
of elevated temperature crack initiation of the hot section structures. 

Ranking Procedure 

A ranking procedure was developed for evaluation and ranking of the life 
prediction approaches and models considered in this task. The ranking 
procedure was based on four criteria considered to be important in assessing 
the usefulness and predictability of a life prediction model. The four 
criteria were: 

A. The amount and types of test data (called baseline data) required to 
determine all of the constants in a model to the same degree of 
confidence. 
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B. The ability of a model to regress the baseline data and accurately 
predict verification test data not used In determination of the 
constants. 

C. The fundamental basis for the model (mechanistic vs phenomenological). 

D. The amount of judgement required to apply a model to an actual 
component loading cycle. 

Initially, a numerical procedure was developed that assigned point values or 
scores for each of the criteria. The first three, A, B and C, were scored 
together, with the first two each given twice as much "weight" as the third. 

The total score for A, B and C was then compared to a score assigned for 
criteria D. This was based on the notion that a model which accurately 
predicted simple cycle data but did not appear to be adaptable to a more 
relevant component loading cycle would not be considered as desirable as a 
model which did not predict the data as well, but potentially captured 
relevant loading cycle characteristics. 

In reality, determination of a consistent score for a number of models 
required that all of the test data be available for each model. This was not 
practical due to the number of tests being run in support of the program and 
tne fact that the test matrix evolved as trends and effects on fatigue life 
were identified. The evaluation process that was eventually used was based on 
the four criteria discussed above, but each model was evaluated with the data 
available at the time of the evaluation. Models that displayed obvious 
limitations with a control data set were not given further consideration, 
while models that displayed certain predictive characteristics were 
continually evaluated as more data became available. In addition, the crack 
initiation life for each test condition changed slightly throughout the course 
of the testing with additional replica data. Ultimately the approach and model 
which developed from the evaluation procedure, and described in detail in Task 
III, represents a combination of attractive features of several of the life 
prediction models considered. 
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3.2.2 Coffin - Manson Model 


The Coffin-Manson Model (Ref. 13) was selected as a starting point for 
evaluation of life prediction approaches. The model ranking criteria, 
discussed in Section 3.3, were used as a "framework" for the initial 
evaluation with available fatigue data grouped into three sets designated A, 
B, and C. 


A; 871 °C (1600°F) , c = 1.67-03 sec" 1 , R = -1 

B; 871 °C (1600°F) , e = 1.67-04 sec' 1 , R = -1 

C; 538°C (1000°F), « = 1.67-03 sec' 1 , R = -1 

Set A was considered baseline data and used for regression of model constants. 
Set B did not show a significant strain rate effect on fatigue life and was 
also considered in the initial model regressions. Set C was used as 
verification data to evaluate the predictive capability of each model. 
Prediction of the lower temperature 538°C (1000°F) fatigue lives with the 
model developed from the higher temperature 871 °C (1600°F) data was ini tally 
considered a valid test of the predictive capability of the model because of 
similarities in the physical initiation site and transgrannular growth mode 
observed at both temperatures. 


The crack initiation life for each specimen test was defined as the number of 
cycles required to produce a 0.76mm(0.030 in.) surface crack. Based on the 
replica data available at the time of the model evaluation, this was 
equivalent to approximately 15 percent of the total specimen life for tests 
conducted at 871 °C (1600°F) and approximately 50 percent of the total life at 
538 # C (1000°F). Previous metallurgical examination had indicated that this 
portion of the life represented a crack depth of approximately two grain 
diameters and that the crack progressed transgranually from the initiation 
site. This observation appeared to be valid for all three data sets. For the 
model evaluation, it was assumed that the crack initiation process was 
represented by the cycle response parameters (inelastic strain and rate, etc.) 
as measured at the crack initiation life. 
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The model assumes a relationship between inelastic strain range and fatigue 
life. Regression and prediction of the three data sets with this model Is 
presented in Figures 42 to 44. 

As shown, combining data set B with set A (Figure 43) did not significantly 
change the model regression. This Is due principally to a small strain rate 
effect on fatigue life for the conditions tested. Using this regressed model 
as a prediction of data set C (Figure 44) resulted In a conservative 
prediction of the crack Initiation life at 538*C (lOOO'F). A later surmnary of 
test data including 538°C (1000°F), 760 # C (1400 # F), 871*C (1600 # F) and 982 # C 
(1800 F) results is shown in Figure 45 and Indicates that a single Independent 
parameter (inelastic strain range) would not correctly present the trends In 
initiation lives. This observation led to the consideration of life prediction 
approaches that Included more than a single Independent parameter. 

3.2.3 Damage Rate Model 

The Damage Rate Model proposed by Majumdar was considered representative of a 
more mechanistic approach and Includes two Independent variables. The model 
assumes that low cycle fatigue Is primarily a process of crack propagation and 
cavity growth. Micro-cracks and cavities are assumed to be present in the 
virgin material and the majority of the low cycle fatigue life Is spent 
growing these micro cracks and cavities to a critical crack size at which time 
they combine to form macro-cracks. 

Following the development presented In Reference 12, the model assumes the 
growth law for the micro cracks and cavities as: 


II 


(For tensile stresses) 

(For compressive stresses) 

(1) 

dc = 
dt = 


(For tensile stresses) 

(For compressive stresses) 

(2) 
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INELASTIC STRAIN RANGE 






INELASTIC STRAIN RANGE. A. in 



Figure 44 Coffln-Manson Prediction of Verification Data (C) 



Figure 45 Variation of Initiation Life with Inelastic Stralnrange and 
Temperature 


where: a = crack length 

c = cavity length 
t = time 

ei - n = inelastic strain 
€j n = inelastic strain rate 
T, C, G, m, K, Kc = constants 

For the current evaluation, a predictive form of the model appropriate to 
continuous cycling at a constant strain rate was used: 

Nf = ( e in) (1 " k) (3 

This form is determined by integrating equations 1 and 2 over an ideal cycle 
which contains equal (tensile and compressive) peak stresses and inelastic 
strains, and constant inelastic strain rate. In the model development, it is 
assumed that the primary cavity growth mechanism is associated with tensile or 
compressive hold times. Thus, the continuous cycling form of the model does 
not include cavity growth. This was considered to be consistent with the 
observed transgrannular process in the tests used for the evaluation. 

Regression of the model (Equation 3) to data set A indicated that the 
two-parameter model was able to better fit the baseline data than the 
Coffin-Manson Model discussed above. A comparison of the baseline data and the 
prediction for various inelastic strain rates (constant total strain rate) is 
shown in Figure 46. 

Application of the model to data set B results in predicted lives that are 
significantly smaller than the experimental data. This is due to the lack of a 
strain rate effect on crack initiation life for the conditions tested (see 
Figure 47). A similar situation exists in the prediction of the 538°C (1000°F) 
tests in data set C (see Figure 48) using the 871 °C (1600*F) model constants. 
The preliminary results obtained with both the Coffin-Manson and Continuous 
Damage Models indicate that the use of single valued macroscopic parameters 
(inelastic strain range and rate) are not sufficient to predict fatigue lives 
at various temperatures. Furthermore, the data used in this initial evaluation 
was developed at zero mean strain (Re = -1). The predictive capability of each 
model was not evaluated on non-zero mean strain. 
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Q 538°C (1000°FI * = 1.67E 03S’ 1 



CRACK INITIATION LIFE 
AT 1 600 °F 0 1 5 * SPEC LIFE 
AT 1000°F 0 60 x SPEC LIFE 


Figure 48 Prediction of Data Set C by Damage Rate Model 


3.2.4 Tensile Hysteresis (Ostergren) Fatigue Model 

General observations of the complete fatigue data set and experience with two 
previous models suggested that a single Independent parameter (Ae1n,Aa, 
etc.) Is not sufficient to correlate the fatigue life at all test conditions. 
It appeared that some combination of stress and inelastic strain would be a 
better parameter. To Investigate this assumption, a data set of nine tests run 
at either 871 °C (1600°F) or 932°C (1800°F) with various strain rates was 
selected. As a baseline, a regression of the data using inelastic strain 
range, as the Independent parameter, and specimen separation life was 
performed and is shown in Figure 49. The data tends to segregate with 

p 

temperature and is reflected In the poor correlation coefficients (R = 

0.77, SEE = 0.23). Considering the product of the stress range and Inelastic 
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□ 1600°F ( 

- 1.67E-3 SEC-1 

871 °C ( 1600°F) 

■ 1600°F * 

* 1.67E-3 SEC-1 

+ 1 MIN COMP. HOLD 

982°C(10OO°F) 

(J 1600°F i 

- 8 3E 5 SEC-1 


y 1 600°F i 

= 1 ,676-4 SEC* 1 


@ 1800®F ( 

- 1 .67E-3 SEC-1 


(Q) 1600 8 F i 

- 1.67E-4 SEC-1 




Figure 49 Life vs Initial Inelastic Strain Range 


strain range (Actx Aein) as the correlating parameter, results in improved 
regression coefficients, as shown in Figure 50 (R^ = 0.90, SEE = 0.15). The 
871 °C (1600°F) is slightly more ordered with respect to strain rate and not as 
much segregation of the 932 °C (1800°F) data is observed. Following an approach 
analogous to the Ostergren Fatigue Model [14], it was assumed that the 
hysteresis loop area might be a better correlating parameter. This results in 
the regression shown in Figure 51 which has the best coefficients of the 
models considered (R^ = 0.93, SEE =0.13). In this correlation, the 
hysteresis loop area was calculated for each test and the cycle frequency (as 
suggested by Ostergren), not included in the model regression. This parameter 
provided the best ordering of the 871 °C (1600°F) data considering three strain 
rates and the effect of a one minute compressive dwell. In addition, the 
smallest amount of temperature segregation was observed. Based on the 
preliminary results, further investigation of this approach to a larger data 
set appeared worthwhile. 
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TOTAL HYSTERISIS LOOP AREA. KSI Aa x AdN. KSI 
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Figure 50 Life vs Stress Range x Inelastic Strain Range 
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Figure 51 Life vs Hysteresis Loop Area 
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Evaluation of the tensile hysteresis loop area as a correlating parameter was 
then evaluated on a larger data set. A baseline data set consisting of rapid 
cycle (e=1.67E-03 S' 1 ), fully reversed (R=-l) tests at 760 # C (1400°), 871°C 
(1600°F), and 932 °C (1800°F) (20 points) was established to determine 
regression constants. Two models were regressed, the first considered only the 
871 °C (1600°F) baseline data (5 points): 

N 1 = 36932 (Area) -1 ’ 504 (4) 

(R 2 = 95.1, SEE = 0.21) 

The second considered the entire baseline data set (20 points): 

N 2 = 14592 (Area)" 1 * 1704 (5) 

(R 2 = 84.5, SEE = 0.31 ) 

Life (N) in equations 4) and 5) is defined as the crack initiation life. The 
"Area" term was calculated as the average of the initial and stabilized 
tensile hysteresis response. A verification data set consisting of 20 data 
points was used to evaluate the predictive capability of both equations. This 
data set was specifically selected to evaluate the uniqueness of the tensile 
hysteresis area for creep fatigue and non-zero mean strain conditions. Test 
conditions included in the data set were: 


1. 871 °C ( 1600 °F) , R e =-1 , e = 1.67E-04, 0.833E-04 S 

2. 871 °C (1500°F), R e =-1, 1 minute compression hold 

3. 871°C ( 1600°F) , R e =-1» 1 minute tension hold 

4. 982°C ( 1800°F) , Re=-1, 1 minute compression hold 

5. 760°C (1400°F) , R e =0, e =1 .67E-04 S" 1 

Correlation of the actual vs. predicted lives for conditions 1, 3, 5 is shown 
in Figure 52. The consistent under-prediction of the lives suggests that the 
total tensile hysteresis area is not damaging. This Is particularly evident 
for the tensile hold (3) and Re=0 (5) conditions where a significant amount of 
mean stress relaxation occurs. A portion of the cyclic response must be 
non-damaging or anelastic based on the above observation. Prediction of the 
complete verification data set by equation 5 is shown In Figure 53. Again the 
general trend is to under-predict the actual lives. 
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2 



ACTUAL LIFE 


Figure 52 Prediction of 871 °C (1600°F) Verification Data Set by Tensile 
Hysteresis Loop Area 



ACTUAL LIFE 


Figure 53 Prediction of Total Verification Data Set by Tensile Hysteresis 
Loop Area 



An estimate of the maximum error In calculating the hysteresis loop area was 
made by repeatedly analyzing a 538°C (lOOO’F), c «+o.25% test. This represented 
the narrowest hysteresis response encountered to date. A calculated variation 
in the area of +20% from the mean value was observed. This was considered to 
be too large a variation for life prediction evaluation. Based on the 
predictive capability discussed above, and the uncertainty In the area 
calculation, this life prediction approach was not considered a viable 
candidate for future consideration. 

3 . 2.5 Local Crack Tip Parameter Models 

The experience with the three models (approaches) discussed above suggested 
that the "far field" macroscopic hypteresls response may not be representative 
of the local conditions which result In creep-fatigue crack Initiation. This 
resulted in the evaluation of two approaches In which local crack tip 
conditions were used as correlation parameters for fatigue Initiation. The 
results reported in [15] Indicated that crack tip opening displacement (CTOD) 
was partially successful In correlating elevated temperature fatigue crack 
growth. As a preliminary assessment for this work, it was assumed that the 
fatigue life could be correlated with CTOD, l.e.: 

N = K (CTOD)^ (6) 

where the McMeeking [16] relationship was used to define the CTOD as: 


CTOD 


0.55J T 2 

[vr 


(1 +»>) (1 + n) 



(7) 


where: n 

O’ 

i 

E 

v 


strain hardening exponent In stress-strain representation 
yield stress 
Young's Modulus 
Poisson's Ratio 
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Calculation of the J integral considered the sum of elastic and plastic 
components following Shih and Hutchinson [17] and assumed that the tensile 
portion of the cyclic hysteresis loop could be used to estimate the cyclic 
parameter. 

( 8 ) 
(9) 


Initial calculations of J at 871 °C (1600*F) for fully reversed (Rt=-1) strain 
range of +0.4%, +0.25% and +0.2%, shown In Table I, Indicated that Jpl is less 
than 30% of the total J. 


J ■ J el ♦ J P 1 




n+1 hJjF "maxi"* 1 


~T~o o"J 


where: h^,rj = functions of crack geometry 

a =1/2 crack length 

a,n = parameters in Ramberg Osgood stress strain 
representation 
°o ■ yield stress 

c = yield strain 


TABLE I 


PREDICTED J AT 871 °C (1600°F) 


At (%) Jel Jpl 


+0.4 

0.0203 

0.0070 

+0.25 

0.0100 

0.0004 

+0.20 

0.0075 

0.0003 


Therefore, for this exercise, J was calculated as: 


J 


» Jel » 



( 10 ) 


67 



For the regression it was assumed that the CTOD associated with the 
"initiated" 0.030-inch crack characterized the entire growth process up to 
this crack length. Regression to a 871°C (1600*F) data set (R€=-l, e=1.67E-03 
S" 1 ) with an "initiation" life associated with a 0. 030-Inch surface crack 
resul ted i n the model : 


N = 5.72E-1 1 (CTOD) -2,94 (11 ) 

Experience with this model is summarized in Figure 54 where the predicted and 
actual initiation lives for the baseline data set and two validation data sets 
(871 °C (1600°F), Re =-l , e=1.67E-04 S' 1 and 538 # C (1000°F), Re=-1, c=1.67E-03 
S _1 ) are shown. As indicated, the predicted lives for both validation data 
sets are consistently longer than the actual data. For the 871 °C (1600°F) slow 
strain rate tests, this is the result of a reduction of the crack tip J (CTOD) 
due to a reduction of the maximum tensile stress relative to the yield stress 
(amax vs. a o in equation 10). This would appear to be contrary to the time 
dependent damage and subsequent reduction on life of the slower strain rate. 
For the 538 # C (1000°F) data set, the predictions also suggest a reduction in 
the predicted value of J; however, this is associated with the small 
difference between the maximum tensile stress and yield stress (o-max vs. <*o in 
Equation 10) consistent with a reduced degree of plasticity at 538°C (1000°F) 
relative to 871 °C (1600°F). Based on these two predictions, it appears that 
the single parameter J is clearly underpredicting the local damage which 
produces initiation. 

Use of a "classical" linear elastic fracture mechanics for prediction of the 
initiation life from a small flaw (carbide or porosity) to the 0.76 mm 
(0.030") surface crack length was also considered. For the two-dimensional 
surface crack the value of the stress intensity factor (K) was determined 
based on finite element results presented in [18]. As a conservative approach 
the value of K at the free surface was assumed to exist along the entire crack 
face and was calculated as: 


AK 


2.4 


a max «7ira 


( 12 ) 
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Figure 54 Correlation/Prediction of Life to a .030 Surface Crack by a 
"Characteristic" CTOD 


which assumes that the crack opens at zero stress. Application of this 
approach together with available crack growth rate data consistently produced 
predicted lives much longer than the actual lives. Average predictions at 
temperatures for fully reversed ( Re=-1 ) rapid cycle testing (Figure 55) 
indicate that better predictions were obtained at the higher temperatures and 
for crack growth from 0.015 to 0.100" In half length, however the predictions 
are at least a factor of 2 longer than the actual average life. A portion of 
the overprediction may be associated with the fact that theAK levels of the 
small crack size was below the range of the test data used to determine the 
growth rate, thus there may be some uncertainty in the actual growth rate at 
low values of AK. It was also judged that generation of reliable crack growth 
rate data for the relevant AK levels associated with the Initiated crack size 
represented a significant burden for development of a life prediction model. 

Based on these results, the use of a CTOD parameter and linear elastic 
fracture mechanics were not considered as practical approaches for the 
prediction of creep-fatigue initiation of B1900+Hf. 
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Figure 55 Best Prediction at Larger Crack Size and Hotter Temperatures 
3.2.6 Total Strainrange Partitioning 

All of the approaches (models) evaluated up to this point have used a single 
correlating parameter. The experience discussed above has indicated that these 
models have not been able to properly account for the relevant effects 
associated with creep-fatigue interaction (strain rate and strain hold time 
effects), mean strain and temperature. Therefore, several models that 
explicitly included separate damage and mean stress correction terms were 
considered. The first of these was based on the concept of strain range 
partitioning. 
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The concept of strain range partitioning [3] assumes the existence of separate 
time independent and time dependent damage life relationships and includes a 
method of predicting fatigue life from a combination of these independent 
relationships. Most of the early application and success of this approach has 
been with materials and loading conditions in which the critical cyclic 
hysteresis loop included a measurable amount of inelastic strain. For cycles 
in which the cyclic inelastic strain is small or non-existent, Halford and 
Sal tsman [19] have proposed a method which extends the concept to the use of 
total strain-range. Following the above method, and including an explicit mean 
stress correction as proposed in [20], the life relationships for a B1900 + Hf 
871 # C (1600°F) data set were determined with subsequent application to a 
separate data set for determination of the predictive capability of the 
approach. The method is based on the decomposition of total mechanical strain 
into elastic and inelastic (plastic) components, i.e.: 

^ elastic +i ‘inelastic 031 

It is assumed that the fatigue life relationships can also be combined to^ 
represent a relationship between total strain range and life. For the 871°C 
(1600°F), Re=-1 , rapid cycle data set, inelastic and elastic life 
relationships for time independent initiation life (pp m strain-range 
partitioning nomenclature) were determined as: 


Ac in = 0.291 (Npp) -0 * 892 


A eel = 0.0203 (Npp) 


-0.179 


(14) 


(15) 


The total strain-range relationship is now written as. 

= 0.0203 (Npp) -0 * 1 79 + 0.291 (Npp)" 0 * 892 


Ae 


TOTAL 


(16) 
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and is shown in Figure 56. In a similar manner, time dependent life 
relationships were determined with the tensile and compressive strain hold 
tests. These tests were conducted with a transient strain rate equal to the 
rapid cycle tests but included a one minute strain hold period. This was 
considered representative of pc and cp fatigue damage in strain-range 

partitioning nomenclature. The life relationships were determined from the 
test data as: 


for pc: 
for cp: 


Ac 7QTAL P c = 0*0185 (Nf ) 
Ac TOTAL c P = 0.0185 (Nf) 


"°* 179 + 0.407 (Nf)" 0,892 

(17) 

"°* 179 + 0.698 (Nf)" 0,892 

(18) 


and are shown in Figure 57. These equations reflect the assumptions that all 
of the inelastic strain observed in the tests was either pc or cp (not pp) In 
nature and that the elastic life relationship was the same for both tensile 
and compressive hold tests. 



Figure 56 Strain Range Partitioning Life Relationships for B1900+Hf Tests 
( e = 1.67E-03 S' 1 ) 
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LIFE. CYCLES 

Figure 57 Time Dependent Life Relationships Used for Total Strain - 
Strainrange Partitioning Model 

The above equations were derived for tests whose duration was 66 sec (6 second 
transient and 60 hold). The elastic portion of the life relationship Is 
considered to be time dependent. An approach used to adjust the relationship 
with time was proposed in [19] and used here. Basically, the elastic life 
relationship is assumed to be written as: 

A eel, ij = Bij (Nij) C H9) 

where ij = refer to the various SRP cycle types pp, pc, cp, cc. 

A time dependent relationship between the coefficients Is assumed as: 

In (Bpp/Bij) * Aij (t) a (20) 

From the available data, constants were determined for the above equation as: 

0.25 (2i) 


In (Bpp/Bij) = .0333t 
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where t = holdtime, seconds. The exponent (0.25) Is based on AF2-IDA 
and Rene '95 data reported in (19). 

The elastic portion of the total life relationship can now be determined for 
various cycle periods. 

Evaluation of the predictive capability of this life approach considered three 
different data sets as verification data. These are: 

1. 1600°F, R € =-l, e =1 . 67E-04 S" 1 

2. 1600°F, R e =-1, €=8.33 -05 S' 1 

3. 1600°F, R € = 0, t=l .67 -03 S" 1 

Prediction of the continuous cycling tests (1 and 2) was made with two 
different assumptions concerning the nature of the damage. In the first 
calculation, the damage was assumed to be equally split between pc and cp 
damage. 

In the second, it was assumed that all of the damage was pp in nature. A 
comparison of the predicted vs. actual lives for both calculations Is shown in 
Figure 58. The second method of calculation (all pp damage) produces a 
slightly better correlation but both sets of results fall within the scatter 
band of the regression data used to determine the model constants. Prediction 
of the Rt=0 tests (3) considered that the damage is characterized as pp in 
nature with equation 16 representing the total strain range life relationship. 
Tnis represents a zero mean stress equation since fully reversed (R«=-l ) data 
was used for the model constants. Adjustment of the life for the observed mean 
stress used the approach discussed in [20] which states that the effect of 
mean stress on fatigue life can be expressed as: 


where: = fatigue life at observed mean stress 

Nf 0 = fatigue life at zero mean stress 
am = observed mean stress 


aa = cyclic stress amplitude 

b = exponent determined from N. = A(Aa) - ^ 

fo 


( 22 ) 
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A comparison of the predicted vs. actual lives indicates a consistent 
underprediction by the model (Figure 58). These Re=0 tests were somewhat of an 
anomaly for most of the life models considered since the data indicates an 
increase in the median life relative to the Re=-1 tests. The approach used for 
the mean stress correction would clearly predict a decrease in life with 
increasing positive mean strains and contributed to the conservative 
prediction of these tests. Despite these predictions, the total strain range 
approach demonstrated encouraging positive capability on a complex 
verification data set. The incorporation of two damage parameters and a 
separate mean strain correction are responsible for the improved predictive 
capability. The disadvantage of the approach is the need for complex baseline 
data to determine model constants. The strain hold tests required for the pc 
and cp coefficients are lengthly (costly) and generation of a sufficient 
number of data points for development of a fatigue design system was 
considered an unattractive data requirement. Therefore, the method was not 
considered as the best overall life prediction approach but several of its 
characteristics were considered in future model evaluations. 


PREDICTED 

LIFE 



Figure 58 Prediction of Continuous Cycling and Re= 0 Tests by Total Strain 
Range Partitioning Method 
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3.2.7 Ductility Exhaustion Model 


One of the less desirable aspects of the total strain-range Partitioning model 
discussed in Section 3.2.6 was the requirement for complex strain hold testing 
to determine the time dependent damage component. This led to consideration of 
models in which the time dependent damage component could be determined from 
monotonic creep data. Models which fell into this category are Ductility 
Exhaustion [21] and Chaboche's Continuous Damage [22, 23, 24]. Fundamental to 
the Ductility Exhaustion model is the assumption that a decrement In available 
ductility, resulting from a period of fatigue cycling, can be represented 
directly as an increment of creep damage. The fatigue component is represented 
as a log-linear relationship between total strain range and life: 


where: 


Ac t - K( Nf )- ,/a 
r * i 1/a 

k *[( 2£ f) /4 ] 


(23) 


2€ F = 1/4 cycle ductility 

The decrease in available ductility produced by one fatigue cycle is then 
written as: 



(24) 


Assuming that creep increment can be described by a power law relationship as: 

ie CR * Att> 25 » 

The equivalent fatigue time increment Is 

*f -W' /b ,2s 

the decrease In the total ductility associated with creep and fatigue damage 
assumes that the time increments can be combined as: 


FR 


= e 


Fi 


A(tf + t creep ) 


(27) 
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where: = residual ductility at end of increment i 

e Fi = initial ductility at beginning of increment i 

Equations 24, 26 and 27 are evaluated on a cycle by cycle basis until the 
original ductility is exhausted. Because of the combined time increments 
raised to a power, the exhaustion process is highly non-linear. 

This approach was originally developed for a situation of gross section (load 
controlled) creep super imposed on fatigue cycling. Application to an 
integrated time dependent (creep- fatigue) strain controlled loading cycle 
required several assumptions on the interpretation of the creep damage 
equation. These cycles represent a variable stress loading cycle. In addition, 
the times associated with the various fatigue cycles are significantly shorter 
than the typical time normally associated with steady state creep. For 
generality, the creep law used in the model was written as 

€ cr = C A 8 *" (28) 

where B, C, m, n are temperature dependent constants. During the course of the 
evaluation, several measures of the stress and time values of the fatigue 
cycles were considered. As an example, consider the prediction of the data set 
shown in Table V below. 


Table V 

Data Set for Evaluation of Ductility Exhaustion Model 
871 °C (1600°F), R € =-1,A€=0.5S 

e = 1.67E-04 sec -1 7 points 

e= 0.83E-04 sec -1 3 points 


From the observed cycle hysteresis loops, the value of the stress for 
determination of the creep damage was calculated as: 


a eff = 


T - 20 


Where: » maximum tensile stress at 1/2 


(29) 
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This represented an average tensile stress for these cycles. A time increment 
equal to 1/2 of the cycle period was used to calculate the creep damage 
Increment per cycle. As shown In Figure 59, the prediction agrees closely with 
the tests run at a strain rate of €=1.67 -04 S" 1 while the slower tests 
(«=.83E-04 Sec~^ ) are consistently overpredicted. These results suggested 
that the increase in time dependent (creep) damage per cycle was not 
sufficiently predicted by equation 28, the average tensile stress and the 
larger cycle period alone. Use of a larger stress magnitude would Increase the 
time dependent damage and improve the correlation between predicted and 
actual lives for these cycles, but would also reduce the predicted lives for 
tne faster strain rate test points as well. Further evaluation of the stress 
for calculation of the creep damage considered the Larson-Miller Parameter. 



Figure 59 Prediction of Continuous Cycling Tests by the Ductility Exhaustion 
Approach 

In Figure 60, the Larson-Miller parameter (P) for all fully reversed fatigue 
tests are plotted vs. maximum tensile stress. The time constant (t) used in P 
represents the time to a 0.030“ crack as determined from replica data. Also 
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MAXIMUM TENSILE STRESS, MPA (KSII 


Figure 60 Larson-Miller Curve Representation of Fatigue Data Using Maximum 
Tensile Stress 

shown in the figure are the creep curves representing 0.2%, 0.4% and rupture 
as determined from the monotonic creep tests. The 0.2% and 0.4% creep limits 
are considered to represent the approximate magnitude of the grain "ductility 
that is exhausted prior to transgrannular crack initiation. As shown, the 
higher temperature tests clearly fall within the creep curves while the lower 
temperature tests have a lower P parameter relative to the creep lives. Not 
surprisingly, this indicates that the fatigue lives are less influenced by 
creep at lower temperatures. In Figure 61, the same data is shown, however, 
the stress associated with each test is the average of the mean and maximum 
tensile stresses. The reduction in the stress magnitudes shifts the points 
away from the creep curves and would be consistent with a much smaller grain 
"ductility". Varying the definition of initiation time (t) does not have a 
significant effect on the position of the data points relative to the creep 
curves as a change in the parameter (P) of 2 points is associated with an 
order of magnitude change in the time. Based on these observations and the 
assumption of the magnitude of the grain ductility, the maximum tensile stress 
appeared to be the more appropriate stress level for calculating the creep 
damage increments. 
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Figure 61 Larson-Miller Representation of Fatigue Data Using Average Tensile 
Stress 

Further evaluation of the model was based on the expanded data set shown in 
Table VI below. This was considered representative of the life trends 
throughout the specimen test program. Relative to the fully reversed, rapid 
cycle tests, life decreases at lower strain rates. Is reduced by a compressive 
strain hold, but Increases with a tensile strain hold. The life also generally 
decreases with increasing (more positive) R ratio with the single exception 
between R = -1 and R = 0. This effect was associated with a longer nucleatlon 
life at R = 0, as a result of the denser dislocation structure produced by the 

loading to a higher value of stress and Is discussed In more detail In Section 
3.3.2. 

Application of the Ductility Exhaustion Model to this data set Is shown In 
Figures 62 and 63. The fatigue model component was determined from fully 
reversed (R * -1 ) rapid cycle (e= 1.67-03 Sec -1 ) data. It was assumed that 
the creep damage component was produced by the maximum tensile stress In the 
cycle acting through one-half the transient period plus any strain hold 
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periods. Application of the model to the prediction of a data set, consisting 
of the slower continuous cycling and strain hold tests from Table VI, is shown 
in Figure 62. The overall trend in the life is predicted although the 
consistent underprediction is most likely associated with a too large creep 
damage component. The model also does not differentiate between the fully 
reversed tests (R = -1) run at two different frequencies (1 cpm vs. 0.5 cpm)*. 

identified by the filled in symbols 


TABLE VI 

SUMMARY OF MEDIAN FATIGUE LIVES AT 871 °C (1600°F) 

Atj = 0.5% 




I 

Life to 

Life to 



Test Condition _ 0. 

030" Crack 

50% Load Drop 

Variable 

R = -1 

1.67-03S" 1 

1250 

3300 

Strain Rate 
1 

R = -1 

1.67-04S- 1 

920 

2750 



R = -1 

8.33-05S- 1 

690 

2050 

1 

( 

R = -1 

1.67-03S- 1 

570 

1700 

Hoi < 

1 Times 


Compression Hold 





R = -1 

1.67-03S- 1 

1290 

3870 




Tension Hold 





R = 0 

1.67-03S- 1 

1780 

4000 

R ratio 

R = 0 

8.33-05S" 1 

510 

1525 

R ratio 





(Strain Rate) 

R = 1/3 

1.67-03S- 1 

1260 

2970 



R = -oo 

1.67-03S- 1 

3040 

6750 

1 

f 


Application to the rapid cycle tests run at three different R ratios 

(0, +1/3, -oo ) from Table VI is shown in Figure 63. The mean strain effect is 

clearly not predicted by the model. The lack of predictive capability was 


expected since the model formulation does not contain a specific "mean strain 
correction term" or recognize the difference in maximum stress and stress 
amplitude of these tests relative to a fully reversed test. Prediction of the 
slower rate (creep dominated) tests but not the rapid cycle mean strain tests 
or vice versa is characteristic of most life models containing two damage 
terms. Experience with the Continuous Damage model on the prediction of 
Hastelloy X creep fatigue tests [25] indicated similar predictive capability. 
Here, this model contains a explicit mean strain (stress) correction, enabling 
prediction of these effects; however, the sensitivity of the results to the 
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PREDICTED LIFE 


871 °C I1600°F) 






stress and time increments for the creep damage made the continued use of the 
model unattractive. These observations led to the model proposed in the next 
section based on the concept of a cyclic damage capability (ductility) that is 
loading history dependent which is exhausted by the independent and time 
dependent damage. 

3.2.8 Summary of Model Evaluations 

The experience with the models, discussed in Sections 3.2.2 to 3.2.7 relative 
to the ranking criteria presented in Section 3.2.1, suggested that an accurate 
and useable model for the prediction of crack initiation under creep-fatigue 
loading conditions should have the following features: 

1. Two damage parameters recognizing both time independent and time dependent 
fatigue mechanisms. As such time (t) should be explicitly included as 
opposed to frequency or strain ratio. 

2. An explicit correction for mean strain or stress. This could also be 
included in a function which is loading history (path) dependent. 

3. Use total mechanical strain range, stress (mean, range, tensile, ...) and 
temperature as the key variables. For the material ( B1 900+Hf ) and relevant 
loading cycles considered here, the macroscopic inelastic strain is small 
and would introduce an additional level of uncertainty in the calculations. 

4. Be formulated with damage ratios as opposed to specific damage increments. 
This would make the prediction relative to a well-defined reference 
fatigue condition. Use of ratios would reduce the potential error 
associated with calculation of specific damage, e.g., cyclic time 
dependent damage using monotonic creep rupture data. 

5. Constants for the model should be completely determinable from simple, 
economic tests, e.g., fully reverse fatigue, tensile, and creep tests. 

All of these considerations were used in the development of the model 
presented in the next section. 
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3.3 Task III - Evaluate Best Candidate Life Prediction Approach 

3.3.1 Basis for Proposed Model 

The discussion presented In Section 3.2 Indicated that none of the models were 
able to sufficiently predict all of the trends in life. It was concluded that 
a more accurate approach could be developed, within the data requirements, by 
selecting certain characteristics of the previous models and interpreting the 
available metal lographic analysis results. This would not only provide a model 
for the specific material considered in the base program (B1900+Hf) but 
potentially be extended to other materials through a qualitative understanding 
of the critical deformation and damage processes. The metal lographic analyses 
conducted on many of the fatigue specimens indicated that the initiation 
process (development of cracks up to 0.030" surface length) is primarily 
transgrannular. This effect is clearly a function of the applied strain rate 
and temperature as shown by the failure map constructed from the creep, 
tensile and fatigue tests (Figure 64). From these observations, the 
fundamental assumption for the development of the model is that the critical 
fatigue damaging processes occur within the grains. Conversely, the critical 
process occurs (or accumulates) within the grain boundaries for those failures 
that are intergrannular. Given these assumptions, a substantial amount of 
examination was conducted to establish a qualitative assessment of the 
deformations observed in various fatigue, tensile and creep tests. This 
information was used to develop a procedure for a quantitative measure of the 
grain (and grain boundary) ability to absorb cyclic damage. The crack 
initiation process is then considered to be an exhaustion of the capability by 
the subsequent fatigue cycling. 

3.3.2 Proposed Interpretation of Material Deformation 

The damage capability is assumed to be related to the grain dislocation 
structure and, as such, the entire loading history must be considered in its 
determination. As an example, fully reversed fatigue tests (Re= -1) of 
B1 900+Hf conducted at relevant strain ranges display little macroscopic 
inelastic strain. Examination of specimens shows a sparse dislocation network 
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within the grains, e.g., Figure 65 A, 871°C (1600°F), Ae = 0.5%, c = 1.67E-03 
sec -1 . For this cycle, the initial straining out to maximum or minimum strain 
(+0.25%) results in local dislocations within the grains. The symmetrical 
nature of the cyclic loading does not create additional dislocations and in 
effect "works" the original structure ultimately producing a crack within the 
grain. For an unsymmetric cycle (Figure 65 B, R = 0) the initial loading out 
to the maximum strain limit produces significantly more macroscopic inelastic 
strain consistent with a denser dislocation structure. If the subsequent 
cycling does not exceed the previous maximum stress (shear), additional 
dislocations are not produced and the fatigue damage accumulates on the 
existing structure. This might be seen in Figures 66 (A, B, C) for a cycle 
between zero and +0.5% strain. During the Initial load cycle (Figure 66 A) a 
dense structure is generated consistent with the macroscopic inelastic strain. 
This structure is considerably different than that produced by a fully 
reversed test of the same strain range; e.g., compare Figure 66 C with Figure 
65 A. Subsequent cycling of another specimen for 70 fatigue cycles (Figure 66 
B) shows an intensification of the dislocations but effectively no additional 
networks. Finally, cycling of a third specimen for 3450 cycles (Figure 66 C) 
shows effectively the same structure. Thus, it was assumed that the 
dislocation network which "absorbs" the fatigue damage Is established during 
the initial loading prior to the cycling. 

Fatigue cycles that include a large component of time dependent damage (e.g. 
strain holds or lower strain rates) displayed a grain dislocation structure 
resembling that observed in a monotonic creep test. As an example, Figure 67 
compares the dislocation structures for a 871 °C (1600°F), Ae = o.5% with a one 
minute compressive strain hold vs a creep test and a slow continuous cycle 
test. These observations suggested that 1) the material grains experience a 
significant amount of deformation during creep and 2) a potential relationship 
exists between the dislocation structure that develops during fatigue cycling 
and that which evolves during the creep deformation. Development of a 
quantitative function for the determination of grain cycle capability from 
creep data further considered the fact that the dislocation structure is a 
function of the applied stress level. Figure 65 A and 66 A showed the increase 
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982°C ( 1 600°F), R 


1 


982 °C ( 1 600°F), R = 0 


Figure 65 Comparison of Dislocation Structures at Two Mean Strains 

(1600°F, Ae = 0.5%. R = 0. e = 1.67E-03S" 1 ) 



1 cycle 



70 cycles 


3450 cycles 


O’ 




Figure 


66 Character of Dislocation Structure Determined by Maximum Stress 






in dislocation structure with stress level associated with the first loading 
cycle of a fatigue test. Interpretation of the monotonic creep response 
assumed that the nonlinear or primary creep was associated with grain 
deformation. Evidence that primary creep strain is dominated by grain damaging 
mechanisms is cited by several authors. Min [26] suggests that the grain 
matrix deformation process at high temperature is a non-linear time dependent 
response which is a good description of the primary creep response. Chambers 
et al [27], in studying the effect of prior creep exposure on fatigue life for 
IN738 concludes that pre-exposure to creep strain levels equivalent to the 
point where transition to minimum creep strain rate occurs (e.g., end of 
primary creep response) is the strain level required to seriously degrade the 
fatigue properties. 

A review of the B1900+Hf creep response at 760°C (1400°F), 871°C (1600°F), and 
932°C {1800°F ) indicated that a linear relationship existed between the 
applied stress level and the amount of primary creep strain. These results are 
presented in Figure 68 and indicates that the amount of primary creep is an 
inverse function of the temperature. This would appear to be consistent with 
the notion of relative grain and grain boundary strengths (equicohesi ve 
temperature) and the observation that lower temperature failures are almost 
exclusively transgrannular. Thus, one of the fundamental assumptions of the 
proposed approach is that the grain cyclic capability can be calculated from 
the maximum stress attained in the initial loading cycle and the amount of 
primary creep strain that could have been developed if the stress level were 
held constant. Since the observed primary strain is bounded at lower stress 
levels, it seems reasonable that higher stress levels (produced by cycles with 
larger positive or negative mean strains) would saturate the grain and then 
produce deformation in the grain boundaries. This assumption was supported by 
two observations. At 871 °C (1600°F), stress levels approaching 95 KSI resulted 
in "slip bursts" in the stress-strain response (see Figure 69). This was 
interpreted as an indication of saturation of the grain and discontinuous 
movement of the grain boundaries. At 760°C ( 1 400°F ) , the boundaries are 
considered to be significantly stronger than the grains. In both the tensile 
and creep tests (see Figure 20b), twins were observed in the grains at 
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Figure 68 Primary Creep Strain is a Linear Function of Stress 
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Figure 69 Observation of Slip Bursts at High Stress Levels 


moderate levels of inelastic strain. This was interpreted as a relief 
mechanism in light of the strong grain boundaries. Thus, the determination of 
the grain cyclic capability required a partitioning of the initial deformation 
to reflect both grain and grain boundary deformation. 

At temperatures < 649°C (1200°F) , creep within the context of the fatigue 
cycles is non-existent. Comparisons of the lower temperature fatigue 
dislocation structures indicated some similarities with the monotonic tensile 
tests (see Figure 70). It was therefore assumed that the residual inelastic 
strain (percent elongation) measured after a tensile test could be considered 
as a measure of the grain cyclic capability. This was considered independent 
of stress level within the fatigue test. A representation of the assumed grain 
cyclic capability vs temperature and stress is presented in Figure 71. 


3.3.3 Grain Damage Equation 

Since a grain based capability has been assumed, it seems reasonable that a 
ductility exhaustion (or damage accumulation) approach be utilized for 
modeling the initiation process. 

A primary assumption for the development is that the fatigue damage processes 
which create transgrannular cracks at a given temperature are similar among 
all transgranularly initiated test conditions, and a single test condition can 
be used to set a reference cyclic damage rate. Any other cyclic damage rate 
can then be calculated by ratioing the reference cyclic damage rate to the new 
condition. 


The proposed functional form of the damage accumulation model is therefore: 


Cyclic Capability 



Cyclic Damage Rate 


Where J [Damage Ratios] 


(30) 
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Figure 70 Low Temperature Fatigue and Tensile Deformation are Similar 




INCREASING 

STRESS 


Figure 71 Assumed Grain Cyclic Fatigue Capability 

The reference damage rate is determined from fully reversed (R e = -1) fatigue 
tests and is illustrated in Figure 72. The maximum tensile stress ( a T ) » 
stress amplitude (A a), and observed crack initiation life (N) for these 
cycles are considered as reference parameters. The grain capability (e p) 
associated with the cycle is determined from <Tj Ref and Figure 71. Thus the 
reference damage rate is the only parameter in the proposed model requiring 
fatigue test data and is calculated as: 

dD 1 _ ?pj (31 

dN jRef "JRef 

This damage rate is obviously a function of the total strain range and can be 
represented as: 

dul - A A, B (3: 

Ref T 

Determination of the damage rate for other fatigue conditions evolved from 
examination of the trends in fatigue life with major cycle parameters. As an 

example, a plot of 871 °C (1600°F) rapid fatigue data, at a strain rate of * = 

0.001667 in/in/sec for various mean strain levels, shown in Figure 73(A) 
suggests that time independent fatigue damage is proportional to maximum 
tensile stress (cr T ) and stress range (Aa). 



_ ep(gR) 

Nr 


= f(Ae) 


Figure 72 Fundamental Form of Reference Cyclic Damage 


1 a °Ind a <r Ao- (33) 

N T 

As shown, a fixed value of stress amplitude combined with Increasing maximum 
tensile stress produces lower Initiation life. From 73 B, this Is associated 
with shifting the hysteresis loop to a more positive mean strain which 
Intuitively results in lower life. Conversely, maintaining a constant maximum 
stress, but decreasing the stress amplitude produces an Increase In life. This 
is possible by decreasing the strain range which again Is consltent with 
longer life. Thus, the total time Independent cycle damage expression Is 
assumed to be related to the reference condition as: 



For fatigue cycles in which time dependent damage Is considered to be 
important, a review of the data suggests a different damage parameter. Shown 
in Figure 74 Is a summary of the 871 °C (1600*F) time dependent and rapid cycle 
baseline (R e- -l ) tests. These results suggest that the damage Is 
proportional to maximum tensile stress (<^ T ) and Inversely proportional to 
stress range (A o') and maybe written as: 

1 1 

1 T “ °Dep a<TJ A^ 05) 
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For these tests, a fixed value of stress amplitude together with increasing 
maximum tensile stress is consistent with lower initiation life. This agrees 
with the time independent damage relationship. However, decreasing stress 
amplitude, at a fixed tensile stress, is consistent with lower life rather 
than higher life. This is a result of the time dependent damage processes and 
suggests that the stress parameters do not solely define the damage. 

Because primary creep strain is assumed to be a good measure of grain 
ductility, it was assumed that a creep power law equation correctly 
represented the relationship between stress, time and damage. Starting with 
the primary creep relationship: 

c p = A a B t C (36) 

where C can be a function of stress ( i . e . , C = Dct E ) and replacing with 
the time dependent stress ratios from equation 35, the time dependent damage 
rate term is written as: 





_ 

- 

B' 

C ’ 1 



dD 

dD 


<tT 

/Ao-Ref \ 


( * ) -i 

(37) 


~3T 


Ref 

a j 

1 Ref 

\ A* ) 


IWef/ J 




Where: t = 1 /2 of the cycle transient period plus hold times. This expression 
assumes that cycles in which the time is equivalent to the time of a reference 
cycle do not contain time dependent damage. This is accomplished by the -1 in 
the expression in the right hand side of the equation. The exponents B‘ and C' 
in equation (37) may not be the same as B and C in equation (36). 

Uetermi nation of the exponents B* and C* is based on data generated in 
monotonic creep tests but with the following two considerations. The exponent 
on stress, B', is determined from strain vs time data in which the origin is 
set at the beginning of the constant load position of the creep test, as is 
tne normal convention. This reflects the sensitivity of the inelasic strain on 
the stress. For the exponent on time, C', the same data is used, however, the 
time origin is shifted back to the beginning of the loading history. This 
results in the elastic ( /E) and inelastic strain being included in the 
regression. 
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This has several effects, 1) inelasticity that occurs in uploading at elevated 
temperatures is included, 2) the measure of time is consistent with that used 
in the baseline reference tests, and 3) the time plays a greater role In 
determination of the damage ratio than if the conventional origin were 
selected. During the course of the model evaluation, use of the exponent C', 
determined from the conventional origin did not predict the life trends as 
well as when C‘ was determined using the shifted time origin. Thus, this 
approach was considered to represent the time dependent damage relative to the 
baseline tests. Combining equations 31 and 37 yields the following damage 
accumulation equation. 




A <x \. /<r T A a Ref 

Ao Ref) ' yr Ref Aa 




dN 


(38) 


where a cycle by cycle Integration is performed to account for evolution of the 
hysteresis loop. 


The entire model can now be written as: 



where : 


€ P 

dD 

“HU 


primary creep strain ductility 



= reference cyclic damage rate 


cry = maximum tensile stress 

Act = stress range 

t = time (total cycle period + hold times )/2 
N = cycle number 

R = reference condition 

B‘ = exponent on stress in primary creep power law 

C 1 = exponent on time in primary creep power law with shifted time 
origin 
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This model has been designated as Cyclic Damage Accumulation (CDA). The use of 
the damage ratios has the advantage of not having to evaluate absolute damage 
increments per cycle. This is considered to be important for cycles that 
include significant amounts of time dependent damage due to the uncertainty 
associated with relating the monotonic creep to the cyclic fatigue process. 

The equation represents a linear damage accumulation approach as the reference 
damage rate is assumed constant through the cycle loading and the damage 
ratios are linearly added. This is considered a good starting point in the 
initial model evaluation. In addition, the available replica data suggested 
that the growth process which defined crack initiation appeared as more of a 
linear process (see Figure 75) as compared to the exponential damage growth 
equations used in the Ductility Exhaustion [21] or Continuous Damage Models 
[ 22 ]. 



Figure 75 Comparison of Surface Crack Measurements and Evaluation of Damage 
Parameters 

3.3.4 Model Predictions 

The CDA model of Equation 39 was used to predict the cyclic crack initiation 
lives of various mean strain, strain rate, strain hold, and block loading 
verification tests as depicted in Figure 76. For the predictions, a 0.030" 
surface crack, which is « 2 grains in depth, was chosen to define initiation. 
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VERIFICATION 

O' 



Figure 76 Types of Cycles Used to Demonstrate Model Predictive Capability 

Application of Equation 39 to the prediction of 81900 + Hf 871 °C (1600°F) 
rapid cycle fatigue tests produced the results shown in Figure 77 A. The 
prediction of various strain R ratios (e min/emax) demonstrates the trend of 
longer life with decreasing mean strain. All test conditions are predicted 
within a factor of +2. For test conditions which are not fully-reversed (i.e., 
R + -1), a ductility debit was determined in evaluating the initial grain 
capability to account for the non-reversed inelastic strain created in the 
first cycle. The method used to determine the ductility debit is discussed in 
Appendix II. 
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77 .lean Strain and Rate Effects Predicted at 871 °C (1600°F) - Model 


Verification Tests 





The ability of the model to distinguish between different loading cycles Is 
further demonstrated In Figure 78. Shown Is the measured crack lengths (from 
replicas) vs. cycles for four (4) types of 871 ’C (1600 # F) rapid cycle tests 
with a strain range of 0.5%. The data suggests that In the fully reversed 
(Ke— 1) test the cracks essentially start from the first load cycle. For the 
other three loading types (Re--, Re«0, Re=0.33), there Is a distinct period 
of little or no crack generation followed by growth at rates that appear to be 
consistent with maximum tensile stress (l.e. Re= .33 has the greatest rate, 
while Re=- * the lowest). Predictions with the proposed model for the 
generation of a .76mm (.030") crack are ordered consistent with the data. The 
net effect of the different time Independent damage ratios and grain 
capabilities, determined for each cycle, may be seen In Table VII where 
calculated values for four representative specimens are presented. As shown, 
the Initial grain capability varies as a function of the maximum absolute 
value of stress generated on the first cycle. The damage ratio of the R — 
test, which reflects the minimum tensile stress. Is less than the fully 
reversed (R =-l ) test due to the compressive mean stress while the ratio for 
the R =0 and 1/3 tests are greater. It is seen that the cycle with the 

greatest grain capability also had the largest damage ratio resulting In the 
lowest 1 ife. 
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TABLE VII 

VARIATION OF DAMAGE RATIOS AND PREDICTED LIFE 
FOR DIFFERENT R RATIO AT 871 °C (1600°F) 


Spec. 

ID. 

1 

Re 

A€ t 2 

e 3 
P 

dD 4 

^Ref 

5 

DRi 

6 

DR 1/2 

N Pred 

23D 

- oo 

.492-02 

.295-02 

.171-05 

0.42 

0.64 

2560 

6 A 

-1 

.50-02 

.192-02 

.185-05 

1.00 

1.00 

1165 

32C 

0 

.497-02 

.353-02 

.179-05 

1.31 

1.04 

1726 

24A 

1/3 

.494-02 

.375-02 

.174-05 

1.69 

1.11 

1141 


2. Total Strain Range, 871 °C (1600°F) , e =1 .67-03 Sec" 

3. Grain Capability 

4. Reference Damage Rate 

5. Time Independent Damage Ratio on First Cycle 

6 . Time Independent Damage Ratio at Half Life 

7. Predicted Initiation Life 


Application of the model to the prediction of 871 °C (1600°F) cycles which 
contain time dependent damage is shown in Figure 77 B. These cycles include 
both rapid transient straining with a one minute tensile or compressive hold 
period and continuously cycling tests with strain rates that are 10 or 20 
times slower than the rapid cycle baseline tests. As shown, the current trends 
in life are predicted. These results can be further examined in Table VIII 
where results for one representative specimen of each test type are presented. 
Relative to a rapid cycle baseline test (specimen 6 A), the model predicts the 
presence of both time independent and time dependent damage. The combination 
of the two ratios represents the total ratio of damage accumulation. 

Initially, all four test types are more damaging than the rapid cycle test. At 
half life, most of the tests still have a combined damage ratio greater than 
one, but the tensile hold test, due to mean stress relaxation, results in a 
damage rate that is smaller than the reference (R=-l) condition. 
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TABLE VIII 

VARIATION OF DAMAGE RATIOS AND PREDICTED LIFE 
FOR TIME DEPENDENT CYCLES AT 871 *C (1600*F) 


Spec. 

ID. 

Ae 

T 

€ 1 

dD 

■3%ef 

Ti/t d i 2 

Ti/T D l/2 3 

( 

Npred 

6A 

.50-02 

1.67-03 

.185-05 

1.0/0 

1.0/0 

1165 

371) 

.50-02 

1.65-04 

.185-05 

.79/. 32 

.75/. 28 

855 

51 A 

.506-02 

833-05 

.196-05 

.76/. 43 

.67/. 44 

728 

27B 5 

.503-02 

1.67-03 

.191 -05 

.89/. 58 

.95/1.11 

479 

22C 6 

.498-02 

1.67-03 

.182-05 

.99/. 69 

.63/. 10 

1455 


1. Strain Rate Sec"^ (A€y Qta ^ /Total Period) 

2. Time Independent Ratio/Time Dependent Ratio on First Cycle 

3. Time Independent Ratio/Time Dependent Ratio at 1/2 life 

4. Predicted Initiation Life 

5. One Minute Compressive Hold 

6. One Minute Tensile Hold 

Extension of the model to the prediction of fatigue tests conducted at 538°C 
(1000°F) and 932°C (1800°F) is presented In Figure 79 (A and B). At 538 # C 
(1000°F) time dependent damage Is not considered and the reduction In life 
associated with positive mean strain Is correctly predicted. The variation In 
life with strain rate at 932°C ( 1 800°F ) Is also predicted. While the 932°C 
(1800°F) results look encouraging, examination of the failure sites Indicated 
an intergranular initiation mode. The form of the model used for these 
predictions is based on the exhaustion of a grain capability and theoretically 
would assume a transgranular mode. Modification of the model to address the 
failure mode is discussed In Appendix IV. 


104 



PREDICTED LIFE 


538°C (1000°F) 



982°C (1800°F) 



Figure 79 Similar Predictive Capability Delayed at Other Temperatures 




6 




To date the model has been applied to simple cycle isothermal fatigue tests 
and has demonstrated the ability to predict various loading effects, e.g., 
mean strain, strain rate, strain holds. Application of the model to the 
prediction of more complex loading cycles considered two types of sequence, or 
block loading, tests. Initially, four tests were conducted at 871 ®C (1600°F), 
at a strain rate of 1.67E-03 sec" 1 , in which the first cycle was applied 
between 0 and +0.5% strain. The mean strain was then adjusted from +0.25% to 
0% and all subsequent cycling occurred between strain limits (+0.28%) which 
produced a maximum tensile stress and stress amplitude consistent with the 
values that would have been obtained in a test with strain limits of + 0.25% 
without the first cycle overload. This sequence was selected to test the 
assumption that the grain cyclic capability Is determined by the maximum 
stress in the loading cycle. Per the model philosophy, the initial overload (0 
to +0.5%) would result in an increased grain capability which, when coupled 
with a damage parameter consistent with the fully reversed parameter (equal 
maximum tensile stress and stress amplitude), would result In a longer 
initiation life. As shown in Figure 80, the initiation lives of the specimens 
are longer than would be expected without the first cycle overload and the 
model successfully predicts the trend in life. 

In the second block loading test, the initial cycling at 0 to +0.5% was 
continued for 600 cylces before making an adjustment in mean strain from 
+0.25% to 0%, while maintaining the strain range at 0.5%. This series of tests 
investigated the rate of ductility exhaustion with variable loading. As shown 
in Figure 80, the model predicted the correct trend In life, however, the 
predictions are consistently conservative relative to the actual data. A 
probable cause of the conservative prediction is that the cyclic damage 
(dL)/dN) in equation 39 is too large during the fully reversed position of the 
cycling. Further investigation of this portion of the model Is required and 
will be considered in future work. 
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(A) 

BLOCK LOADING 




Figure 80 Prediction of Verification Test Cycles 


3.3.5 Model Sensitivity Study 

An important aspect of any life modelling technique Is an understanding of how 
the prediction capability Is influenced by assumptions and specific variables. 
During development, the following parametric studies were identified which 
provided insight to overall model performance. 

1) Crack size - How does the model prediction capability change for 
various ratios of initiation crack size to grain size? 
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2) Crack initiation definition - How does the prediction capability vary 
as a function of the definition of initiation vs. total specimen 
separation life? 

3) Reference condition - Can the model use an alternative reference 
condition to predict crack initiation lives? 

4) Analysis requirements -How much prediction accuracy is lost If the 
cycle dependent hysteretic response is unavailable? 

To present the results, the Weibull slope and median value of the ratio 
predicted life/actual life (Np/Na) were calculated for the 871 °C (1600°F) data 
set presented in Figure 77 assuming a variation in each of the four 
sensitivity areas discussed above. As shown in Figure 81, a larger slope (/3) 
is Indicative of less scatter in the predictions. The median life ratio was 
included to identify consistent trends for under or over prediction. 

Crack Size 


In the experimental program, special efforts were taken during specimen 
casting to control grain size to insure specimen isotropy. In actual 
application, however, grain size may vary from .005" - .150" depending on the 
various casting constraints. It is, therefore, important to establish how the 
model predicts cracking lives over a wider range of crack to grain sizes. 
Utilizing replica data, life fractions for .015", .030" and .060" surface 
cracks (i.e. crack depths of approximately 1 to 4 grains) were defined to 

stucfy prediction accuracy as a function of crack size. 

The model predictions for a 871 °C (1600°F) data set that included both mean 
strain and creep fatigue conditions are shown in Figure 82. Crack depth to 
grain size in excess of 2 showed little variation In the predictive capability 
of the model. For a crack depth on the order of the grain size, a decrease in 
the Weibull slope and an increase in the median predicted to actual life ratio 
indicate a loss in prediction accuracy. This was attributed to local 
non-isotropic single grain effects. A qualitative estimate of model 
requirements to increase the prediction accuracy for a grain size crack 
indicated that the primary creep ductility would necessarily have to be larger 
and more sensitive to applied stress level. 
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CUMULATIVE 


O" = MEAN STRAIN TESTS 
■£,"= CREEP - DWELL TESTS 


N PRED /N ACTUAL 


Figure 81 Parameters Used in Sensitivity Analysis 



CRACK DEPTH/GRAIN SIZE 


Figure 82 Variation of Predictive Capability with Crack/Grain Size Ratio 


Crack Initiation Definition 


The crack initiation life for each test specimen was determined from replica 
data, l.e. , several specimens of each test type were Inspected to determine an 
average life fraction (cycles to Initiation/cycles to 50% load drop) which was 
then applied to all specimens tested at that condition. The observed life 
fractions were simlliar but not equal for all test conditions. 

To access the change in model prediction accuracy, the 871°C (1600°F) data set 
(i.e. data of Figure 77 (A and B) was recalculated assuming a .030" surface 
crack initiated at a constant 38% of the 50% load drop life for all test 
conditions. This resulted In effectively no change In the Welbull slope and 
median life ratio. This would suggest that specific life fractions are not 
necessary for accurate predictions with the model, however, these results may 
be influenced by the types of loading cycles and the range of life data 
generated in the program. It is recommended that Inspections for definition of 
crack initiation be continued until the sensitivity of the model can be 
further investigated. 

Reference Condi ti on 


The current model assumed fully reversed testing (R = -1 ) as the reference 
condition. If this is not available, an alternative reference condition would 
be required. To quantify the change In the predictive capability with the use 
of an alternative reference condition, the 871 °C (1600°F), R = 0, e = 

.001667 in. /in. /sec test condition was chosen as a new reference and 
predictions made of the remaining 871 °C (1600°F) data. As shown in Figure 83, 
the change in reference condition results In a significant loss of predictive 
capability. A review of the specific test types Included In the data set 
indicated that a large part of the loss of predictive capability was 
associated with the slow strain rate, creep fatigue tests. This trend is 
associated with the stress ratio In the time dependent damage term of Equation 
39. Since the new cr TRef is high, the stress quantity is « 1 for the slow 
cyclic rates. In reviewing Figures 74 and 75, it Is evident that the monotonic 
data responses which allowed the simple ratios of Equations 33 and 35 to be 
written are not present when the new reference condition Is used. Therefore, 
it was considered necessary to use R - -1, fast cyclic rate data as the 
reference condition for the current form of the proposed model. 
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Rf = - 1 Rf = 0 

BASELINE BASELINE 


Figure 83 Variation of Predictive Capability with Reference Condition 


Analysis Requirements 

Within the context of the structural analysis of an actual component, the 
determination of the complete cyclic hysteresis history may be impractical. 
More likely, only the cyclic response for the first one or two cycles and an 
estimate of the stabilized response may be available. 

To evaluate the effect of the cyclic response definition on prediction 
accuracy, the first cycle response and a typical 1/2 life cycle response were 
individually analyzed for the data of Figure 77 (A and B). As presented in 
Figures 84, some prediction accuracy is lost by ignoring the cyclic 
relaxation, however, the model still provides an accurate representation of 
the data trends. 
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Figure 84 Variation in Prediction Accuracy with Cyclic Response Definition 



4.0 Future Work 


In the option portion of this contract the evaluation of creep fatigue life 
models will be extended to more complex loading cycles. Testing, 
metal lographic examinations and model evaluations will include: 

1. Thermomechanical loading cycles 

2 . Multi -axial stress and strain states 

3. Cumulative damage loading cycles 

4. Environmental effects on fatigue life 

5. The interaction of oxidation resistant coatings and the underlying base 
metal on initiation life. 

In addition, life prediction approaches evaluated for the cast material 
(B1900+Hf ) will be applied to the prediction of crack initiation of a wrought 
material. This will give an indication of the degree of "transportibil i ty" of 
the models from one material system to another in which the deformation and 
damage processes differ. A wrought model base alloy (IN718) will be used for 
the alternate material evaluation. 

Trie COA approach described in Section 3 will form the basis for these 
evaluations but other model approaches will be included to address specific 
aspects of the tests. 
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APPENDIX I 

SUMMARY OF B1900+Hf FATIGUE TESTS IN BASE PROGRAM 
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APPENDIX II 

EFFECT OF UNREVERSED INELASTIC STRAIN ON NET CYCLIC CAPABILITY 

The cyclic daiiiage accumulation model of Equation 39 is referenced to fully 
reversed, fast strain rate data (R-l, e- .001667 in. /in. /sec) which 
effectively accounts for the cyclic fatigue damage process. Alternate mean 
strains (R^-l), produce large, relative to the cyclic inelastic strain, 
nonreversed inelastic strains during the initial loading, as shown in Figure 
II-l. It was assumed that this nonreversed inelastic strain should be deducted 
from the initial material cyclic capability value. 

As defined, is a measure of the unreversed inelastic strain. At lower 
temperatures where time dependent deformation (creep) is not active, the 
inelastic strain is considered to occur exclusively within the grains. At 
elevated temperatures, the observed inelastic strain is considered a mixture 
of grain and grain boundary deformation and should be partitioned to determine 
the portion affecting the grain capability. In the initial evaluation, the 
partitioning was based on the ratio of primary and secondary creep rates. From 
creep data, the grain componentf primary creep) is represented as: 


Aa B t C 


where C = D 

and the grain boundary component (secondary creep) is: 


€ 

S 


= aV 


t 


(II-l) 


(II -2 ) 


Now, assume e in of Fig II-l can be represented by primary and secondary creep 
strains. 


e . = e + e c 

in p S 


and from Equations II-l and 1 1 -2 


( 1 1 — 3 ) 
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( 1 1 -4 ) 


. B t Da 
A <7 t 


~r a I 

Aa 8 t 


or 



where a= <r I<1AX of the first cycle 

t = loading time + strain hold times which increase e in 


( 1 1 - 5 ) 


Mow let X = e p /e s 

then e = e + e /l\ 
in p p \ x / 



(II -6 ) 


( 1 1 - 7 ) 


The primary creep strain can then be subtracted from the total grain 
capability (e p ) to estimate a net cyclic capability for calculation of the 

initiation life. 

For tests run at Re ■ 1/3 (.25% to .75%) or - « (-.50% to 0), certain 
limitations of the experimental setup resulted in stress relaxation during the 
initial quarter cycle loading. This is seen in Figure 1 1-2 for a Re = 1/3 test 
in which the relaxation occurs at a strain of 0.25%. For the R e = test 
(Figure II-3), the relaxation is integrated throughout the initial loading. 

The additional unreversed inelastic strain produced in these cycles was also 
considered as a debit to the grain cyclic capability and was therefore 
included in the determination of the net value prior to the calculation of the 

initiation life. 
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Figure II-2 Additional Inelastic Strain Produced in R = 1/3 Cycle 




APPENDIX III 

ALTERNATE FORMS OF THE MODEL 


The complete creep fatigue equation proposed in section 3.3.3 is: 



(III-l) 

This represents an approach based on a knowledge of the strain range and the 
local stress response and is considered a creep-fatigue interaction model. By 
considering two special cases, it can be seen that the model "brackets" 
generally accepted life approaches for both low and high temperature. For 
example, by considering only time independent fatigue and assuming that the 
stress response does not change throughout the cycling, the life equation can 
be written as: 


C = a (Ae) b 




N 


(II 1—2 ) 


where: C = e p represents a constant grain capability 

dul = a(A€ ) b 

d,J J Ref 

Ac= total or inelastic strain range 
a, b = regression constants 

This form is similar to life models proposed for low temperature (time a wave 
shape independent) fatigue e.g. Coffin-Manson. The term containing the stress 
ratio can be considered as a mean stress correction. 


For cycles in which time dependent damage dominants (slow, high temperature 
tests), a limiting value of the equation can be determined assuming: 

1. time independent damage is neglible. 

2. the stress ratio in the time dependent term is approximately 1. 
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3. the cycle time is »the reference time. 

With this assumptions equation III-l can be written as 


C -f* M b (^ e jf * (I11 - 31 

where: Aa= EAc 

t = loading time fo reference cycle 

t = loading time of reference cycle 

this fora! the equation appears as an integration of power law creep damage 
which has been shown to be an accurate life prediction method for very high 
temperature tests. Thus, equation III-l represents a transition between low 
temperature time independent fatigue and high temperature creep dominated 

processes. 
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APPENDIX IV 


One of the stated desirable features of a reliable modeling procedure Is to 
distinguish between transgranular and intergranular crack initiation modes. 

Tne bulk of development and evaluation was performed on the transgranular 
model since most of the test specimens failed transgranularly. However, 

transition towards intergranular initiation was observed at 932 °C (1800°F) for 
slow cyclic rates. 

Prel iminarily , two intergranular models were considered. The first model Is 
similar to the transgranular model given in Equation 39, except that all terms 
and constants are related to secondary creep which was considered to represent 
the grain boundary damaging process. This is referred to as method 1. 



wli6r6 * 

€ S = secondary creep strain ductility as shown in Figure IV-1 

Tfir = reference cyclic damage rate 
JRef 

<r T = maximum tensile stress 

A a = stress amplitude 

t = time 

t = 1/2 (Ramp Time) + tensile hold 
N = cycle number 
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Ref 


reference condition 


13' = exponent on stress in secondary creep power law 

C' ■ exponent on time in secondary creep power law 

The second model considered was simply a secondary creep power law ductility 
exhaustion model. This is referred to as method 2. 

f A<7, 8 t dW = 0 

u - J 7 

where 

e = secondary creep strain ductility 
c s 

<jj = maximum tensile stress 

t * time . 

t = 1/2 (Ramp Time) + tensile hold 

A B h coefficients of e * Atr B (assuming linear secondary creep 
’ s strain rate.) 

N = cycle number 
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Application of these two approaches to fatigue cycles in which intergranular 
initiation was observed has shown the ability to predict the current trend in 
the failure mode, although the exhaustion of creep ductility (method 2) 
appears to give more accurate predictions at this point. These results are 
summarized in Figure IV-1 where the predicted ductility (cycle capability) 
exhaustion for four different cycles are presented. In each figure, the cyclic 
capability vs fatigue cycles are predicted using the grain (transgrannular) 
model (equation 39) and the two intergranular models. In Figure IV-1 (A and B) 
the predicted exhaustion for two 371°C (1600°F) tests are presented. Both 
tests were at a strain range of 0.5%, however, the second test (IV-1 B) was 
run at a strain rate one half of the first test. Examination of the fracture 
surfaces of the specimens indicated a greater tendency toward intergranular 
initiation at the slower rate. The predictions for both cycles show that the 
grain capability is exhausted before the grain boundary capability, however, 
the rate of exhaustion for the slower cycle indicated that a further reduction 
in strain rate would have resulted in a predicted Intergranular initiation. 

The second observation on these two cycles is that method 1 (grain boundary 
form of CDA damage equation) is exhausting the available ductility at a faster 
rate than in the second method. This may have some physical significance as 
the second method would seem to be more appropriate at temperatures dominated 
by creep damage processes. This is seen In Figure IV-1 (C and 0) for two 932*C 
(180Q°F) cycles run at the same strain rate but for two different strain 
ranges. In IV-1 C transgranular initiation is predicted but the capability of 
the grain boundary is being exhausted at a rapid rate by the second method. 
Examination indicated this test to produce transgranular initiation. At a 
larger strain range (Figure IV-D), intergranular Initiation was predicted 
consistent with the observed failure mode. The linear creep exhaustion model 
(method 2) appears to be more representative of the failure process at these 
conditions, than the creep fatigue interaction equation represented by method 1. 

Further testing and model evaluations are required before any final 
conclusions are drawn, but both methods may be useful In distinguishing 
failure mode for various fatigue cycles. 
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Figure IV- 1 Exhaustion of Grain and Grain Boundary Ductilities 


3 



DISTRIBUTION LIST 


Aerojet Liquid Rocket 
Attn: V. Frick 
Mgr, Matl 's Eng 
box 13222 

Sacramento, CA 95813 
AFOSk/NA 

Attn: Maj . David Glasgow 

Mgr, Aerospace Sci 
boiling AF8 
Washington, DC 20332 

Wright Patterson Air Force base 

Attn: Ted Nicholas 

AFWAL/MLLW 

Wright Patterson 

Air Force Base, OH 45433 

Wright Patterson Air Force Base 
Attn: Library 

AFWAL/MLLAM 
Wright Patterson 

Air Force Base, OH 45433 


Arn\y Applied Tech Lab 
Attn: Jan Lane 

DAVDL-ATL-ATP 
Fort Eustis VA 23604 

Array Applied Tech Lab 
Attn: Library 

Fort Eustis, VA 23604 

AVCO - Space Systems Div. 
Attn: Library 

Lowell Industrial Pk. 
towel 1 , HA 01 851 

AVCO - Lycoming Div. 

Attn: Jerry Walters 

550 S. Main St. 

Stratford, CT 06497 

Battel le Columbus Lab 
Attn: Brian Leis 

506 King Ave. + 

Columous, OH 43201 


AiReasearcn iifg Co. 
Attn: Richard Graves 

2526 W. 190th St. 
Torrance, CA 90509 


Boeing 

Attn: Library 

3801 S. Oliver 
Wichita, KS 67210 


Univ. of Alabama 
Attn: Ur. A. E. Carden 

AME Dept. 

Box 2908 

University, AL 35486 

Allison Gas Turbine Oper. 
Attn: Mehmet Doner W8 

P. 0. Box 420 
Indianapolis, IN 46206 

Allison Gas Turbine Oper. 
Attn: Frank Walters 

P. 0. Box 894 
Indianopolis, IN 46206 
Allison Gas Turbine Oper. 
Attn: Donald Vaccari T10 

P. 0. Box 420 
Indianopolis, IN 46206 

Univ. of Arizona 

Attn: Dr. Paul H. Wirsching 

Dept Aero & Mech Eng. 

Tucson, AZ 85721 


Case Western Reserve Univ. 
Attn: Prof. S. S. Hanson 
619 Glennan Bldg. 

10900 Euclid Ave. 
Cleveland, OH 44106 

Univ. of Cincinnati 

Attn: Ur. Donald Stouffer 

Rm. b33 Rhodes Hall 
Dept. Eng. Sci. 

Cincinnati, OH 45221 

Univ. of Connecticut 
Attn: Dr. Eric Jordan 

Mech. Eng. Dept. U-139 
Storrs, CT 06268 

Curtiss Wright 
Attn: Jerome Mogul 

Dir - Matl 's Eng. 

1 Rotary Dr. 

Wood-Ridge, NJ 07075 


DISTRIBUTION LIST (Continued) 


Lmi v. of Layton 

Attn: Dr. Joseph Gallagher 

Research Institute 

Kin. 563 Kettering Bldg. 

Dayton, OH ^469 

Defense Documentation Ctr. 
Canieron Station 
5U10 Duke St. 

Alexandria, VA 22314 


Gen. tlec. to. - ALBG 
Attn: Joe McKenzie 

Technical Programs 
mail Drop H-9 
Evandale, OH 45215 

Gen. Elec. Co. 

Attn: Donald Mowbray 

Bldg. 55-219 
1 River Rd. 

Schenectady, NY 12345 


DeLaval Turbine 
Attn: Harry Gayley 

Nottingham Way 
Trenton, NJ 08602 


Gen. Llec. Co. - AEBG 
Attn: Dr. Mel Roberts 

Mail Drop K-69 
Evandale, OH 45215 


Failure Analysis Assoc. 
Attn: Dr. Clifford Wells 

2225 E. Bayshore Rd. 

Palo Alto, CA 94303 


Gen. Elec. Co. - R&D Ctr. 
Attn: Dr. David Woodford 

Bldg. K-l , Rm 231 M 
Schenectady, NY 12301 


Garrett Turbine Eng. Co. 
Attn : Joe Adams { S03-4Z ) 

111 S. 34th St. 

Box 5217 

Phoenix, AZ 85010 

Garrett Turbine Eng. Co. 
Attn: Lee Hatsch 

111 S. 34th St. 

Box 5217 

Phoenix, AZ 85010 

Gen. Elec. Co. - AEBG 
Attn: Dr. Len Beitch K-221 

Mgr, EM & LM 
Evandale, OH 45215 

Gen. Elec. Co. - AEBG 
Attn: J. H. Laflen 

Mail Drop G-oO 
Evendale, OH 45215-6301 

Gen • Elec. Co. 

Attn: Library 

1000 Western Ave. 

Lynn, MA 01905 

Gen. Elec. Co. 

Attn: Library 

Box 8 

Schenectady, NY 12301 


Gen. Elec. Co. - AEBG 
Attn: Kenriard Wright 

Evendale, OH 45215-6301 

Gen. Elec. Co. - AEBG 
Attn: T. Cook 

Evendale, OH 45215-6301 

Georgia Inst, of Tech. 

Attn: Prof. S. Antolovich 

Fracture A Fatigue Lab 
Bunger-Henry Bldg. 

Atlanta, GA 30332 

Hamilton Standard 
Attn: Gerald Mol ter 

Hail Stop 1-1-6 
Bradley Field Rd. 

Windsor Locks, CT 06096 

Hiboitt, Karlsson A Sorensen, Inc. 
Attn: Dr. Kevin Walker 

35 South Angel 1 St. 

Providence, RI 02906 

Univ. of Illinois 
Attn: Prof. 0. Socie 

Dept. Hech tngrng 
1206 W. Green St 
Urbana, IL 61801 


139 


DISTRIBUTION LIST (Continued) 


I. I. T. Research Inst. 
Attn: Humphries 

10 W. 36th St. 

Chicago, IL 60616 

int'l Harvester Co. - Solar 
Attn: Library 

3200 Pacific Hwy. 

San Diego, CA 92101 

Lockheed California 
Attn: Library 

Burbank, CA 91503 

Lockheed Huntsville 
Attn: W. H. Armstrong 

Mgr, Struct <i Hech 
Huntsville, aL 35806 

Mass. Institute of Tech . 
Attn: Prof. Regi Pelloux 

Rrn 8-2S7 
77 l«iass Ave. 

Cambridge, MA 02139 

McDonnell Douglas 
Attn: Library 

Missiles & Space Div. 

5301 Bolsa Ave. 

Huntington Beach, CA 92647 

Michigan State Univ. 

Attn: Dr. John Martin 

MMM Dept. 

330 Engineering Bldg. 

E. Lansing, III 48824 

NASA-AFWAL Tech Liaison 
Attn: Everett Bailey 

AFWAL/DO 

WPAFB, OH 45433 

NASA Headquarters 
Attn: Dr. M. Greenfield 

RTM-6, ms Div. 

Washington, DC 20546 

NASA Langley Research Ctr. 

Attn: Library 

MS-185 

Hampton, VA 23665 


NASA - LeRC 
Attn: Dr. Robert Bill 

Head, Fatigue Res. Sect 
MS 49-7 

Cleveland, OH 44135 
NASA - LeRC 

Attn: D. E. Sokol oski 

Mail Stop 49-1 
Cleveland, OH 44135 

NASA - LeRC 
Attn: T. J. Miller 

Mail Stop 49-3 
Cleveland, OH 44135 

NASA - LeRC 

Attn: HOST Proj. Off. 

Mail Stop 49-7 
Cleveland, OH 44135 

NASA - LeRC 
Attn: R. H. Johns 

Mail Stop 49-6 
Cleveland, Oh 44135 

NASA - LeRC 

Attn: Or. S. R. Levine 

Mail Stop 105-1 
Cleveland, OH 44135 

NASA - LeRC 

Attn: D. B. Ercegovic 

Mail Stop 6-11 
Cleveland, OH 44135 

NASA - LeRC 

Attn: 0. R. Englund 

Mail Stop 77-1 
Cleveland, OH 44135 

NASA - LeRC 
Attn: J. E. Rohde 

Hail Stop 6-10 
Cleveland, OH 44135 

NASA - LeRC 
Attn: M. T. Saunders 

Mail Stop 3-3 
Cleveland, OH 44135 


140 


DISTRIBUTION LIST (Continued) 


NASA - Lt-KC 

Attn: J. A. Ziemianski 

Mail Stop 40-6 
Cleveland, OH 44135 

NASA - LeRC 
Attn: A. Kaufman 

Mail Stop 49-7 
Cleveland, OH 44135 

NASA - LeRC 

Attn: Dr. R. L. Thompson 

Mail Stop 49-6 
Cleveland, OH 44135 

NASA - LeRC 
Attn: H. A. McGaw 

Mail Stop 49-7 
Cleveland, OH 44135 

NASA - LeRC 
Attn: T. W. Orange 

hail Stop 49-6 
Cleveland, Oh 44135 

NASA - LeRC 
Attn: Dr. L. Berke 

Mail Stop 49-6 
Cleveland, OH 44135 

NASA - LeRC 

Attn: Dr. D. A. Robinson 

Mail Stop 49-6 
Cleveland, OH 44135 

NASA - LeRC 
Attn: Dr. R. Ellis 

Mail Stop 49-6 
Cleveland, OH 44135 

NaSa - LeRC 

Attn: Dr. R. L. Dreshfield 

Mail Stop 49-1 
Cleveland, OH 44135 

NASA - LeRC 
Attri : S. J. Grisaffe 

Mail Stop 49-1 
Cleveland, OH 44135 

NASA - LeRC 

Attn: Dr. R. V. Miner 

Mail Stop 49-3 
Cleveland, OH 44135 


naSA - LeRC 
Attn: M. V. Hatha! 

nail Stop 49-1 
Cleveland, 01! 44135 

NASA - LeRC 
Attn: B. McKay 

hail Stop 49-3 
Cleveland, OH 44135 

NASA - LekC 
Attn: Dr. J. Gay da 

hail Stop 49-3 
Cleveland, OH 44135 

NASA - LeRC 

Attn: Dr. Gary Halford 

MS 49-7 

Cleveland, OH 44135 
hASa - LeRC 

Attn: Marv Hirschberg 

Chief, Fatigue & Fracture 
MS 49-6 

Cleveland, OH 44135 

NASA - LeRC 

Attn: Report Control 

hS S-o 

Cleveland, OH 44135 

NASA - LeRC 
Attn: Library 

MS 60-3 

Cleveland, OH 44136 

NASA - LeRC 
Attn: SaT Section 

MS 5Cl -1 1 

Cleveland, OH 44135 

NASA - LeRC 
Attn: ARATL Office 

MS 77-5 

Cleveland, OH 44135 
NASA - LeRC 

Attn: AFSC Liaison Office 

MS 501-3 

Cleveland, OH 44135 


141 


DISTRIBUTION LIST (Continued) 


NASA Marshall Space Ctr. 
Attn: Library 

Huntsville, AL 35812 

NASA Marshal 1 Space Ctr. 
Attn: Carmel o bianca 

EP <+2, bldg oblo 
Huntsville, AL 35812 

NASA Marshall Space Ctr. 
Attn: Larry Salter 

LP46 , Bldg 4610 
Huntsville, AL 35812 

NASA Marshall Space Ctr. 
Attn: Bryan McPherson 

Huntsville, AL 35812 

NASA Marshall Space Ctr. 
Attn: B Bhat 

Huntsville AL 35812 

NASA Marshall Space Ctr. 
Attn: P. Munafo 

EH22, Bldg 4612 
Huntsville, AL 35812 


Naval Ship R&D Ctr 
Attn: Library 

Code 522.1 
Bethesaa, Mb 20084 

Northwestern Univ. 

Attn: Dr. Julia Weertman 

Dept. Matl's Sci 
Evanston, IL 60201 

0 R N L Box X 

Attn: Dr. Charles Brinkman 

Metals & Ceramics Div. 

Oak Ridge, TN 37330 

O'Donnell a Assoc. 

Attn; Dr. W. O'Donnell 
241 Curry Hollow Rd. 
Pittsburgh, PA 15232 


Pennsylvania State Univ. 
Att: ur. Sam Zamrik 

121 Hammond bldg. 

University Park, PA 16302 

Rensselaer Polytech Inst. 
Attn: Prof. Norman Stolofr 

Dept, of Matl 's Lng. 

Troy, NY 12181 

kockwell Int'l 

Att: Dr. Mike Mitchell 

Rockwell Science Ctr. 

P. 0. Box /I 085 
Thousand Oaks, CA 91360 

Rocketdyne 

Attn: Dr. Bob Jewett AC10 

6633 Canoga Ave. 

Canoga Park, CA 91304 

kocketdyne 

Attn: Ed Larson AClu 

6633 Canoga Ave. 

Canoga Pk, CA 91304 

kocketdyne 

Attn: Jack Lewis AC 10 

6633 Canoga Ave. 

Canoga Pk, CA 91304 

Rockwell Int'l 
Attn: Library 

6633 Canoga Ave. 

Canoga Pk, CA 91304 

Rolls Royce, Inc. 

Attn: Dr. Ashok Thakker 

1895 Pfioenix Blvd. 

Atlanta, GA 30349 

Saridia Labs 
Attn: Jack Hyzak 

Division 8314 
box 969 

Livermore, CA 94550 

Sandia Labs 
Attn: Library 

Albuquerque, NM 87115 


NASA Sci. & Tech. Info 
Attn: Accessioning Dept. 

Box 8757 

Balt/Washington International Airport 
Maryland 21240 


142 



DISTRIBUTION LIST (Continued) 


Southern Research Inst. 

Attn: Library 

2000 9th Ave. S. 

Birmingham, AL 35205 

Southwest Research 
Attn: Dr. Gerald Leverant 

P. 0. Drawer 28510 
San Antonio, TX 78284 

Southwest Research 
Attn: Dr. U. S. Lindholm 

P. 0. Drawer 285 
San Antonio, TX 78284 

Southwest Research 
Attn: Library 

8500 Culebra Rd. 

San Antonio, TX 78284 

Stanford Univ. 

Attn: Prof. Alan K. Miller 

Dept. Matl's Sci & Eng 
Stanford, CA 94305 

Syracuse Univ. 

Attn: Dr. H. W. Liu 

409 Link Hall 
Syracuse, NY 13210 

T R W Inc. 

Attn: Dr. C. Kortovich 

T/M-3367 

23555 Euclid Ave. 

Cleveland, OH 44117 

Teledyne CAE 
Attn: Hugh Gaylord 

Mgr., Explor. Dev. Appl. 
1330 Laskey Rd. 

Toledo, OH 43612 

Teledyne CAE 

Attn: Bryon L. Lewis 

1330 Laskey Rd. 

Toledo, OH 43612 

Teledyne CAE 
Attn: Tom Moyer 

Box 6971 

Toledo, OH 43612 


Teledyne CAE 
Attn: Jerry Walcher 

Box 6971 

Toledo, OH 43612 

Univ. of Tennessee 
Attn: Dr. V. Smith 

Space Institute 
Tullahoma, TN 37388 

V. P. I. & State Univ. 

Attn: Dr. Norm Dowling 

Eng. Sci. & l-iech. Dept. 

Blacksburg, V A 24061 

Williams Research 
Attn: Peter Nagy 

2280 W. Maple Rd. 

Walled Lake, MI 48088 

Westinghouse - R&D Ctr. 

Attn: Robert Johnson 

Bldg. 401, Rm. 2X9G 
Pittsburgh, PA 15235 

Westinghouse - ARD 
Attn: Alfred Snow 

Box 1 58 

Madison, PA 15663 

NASA Headquarters 
Attn: Charles Bersch 

RTM-6 , MAS Div. 

Washington, DC 20546 

N. A. P. C. 

Attn: A. Martino 

P.0. Box 7176 
Mgr. R&T Div. 

Trenton, NJ 08628 

Wright Patterson Air Force Base 
Attn: Richard Hill 
AFWAL/POTC 
AREA B, Bldg 18 

Wright Patterson Air Force Base 

OH 45433 

Ai Research Mfg. Co. 

Attn: Dr. Tekal Nath 
Mail T-42, Dept 93-3 
2525 W 190th St 
Torrance. CA 90509 


143 


DISTRIBUTION LIST (Continued) 


University of Akron 
Attn: Prof. Padovan 

Dept. Civil Eng. 
Akron, OH 44325 

Argonne Nat. Lab. 
Attn: S. Majumdar 
9700 S. Cass Ave. 
Argonne, IL 60439 


AVCO - Lycoming Div 
Attn: Louis Fiedler 
Chief, Natl 's Eng & Dev 
550 S Main St 
Stratford, CT 06497 

Case Western Reserve Univ 
Attn: Prof. A. Chudnovskl 
205 Bingham Bldg. 

10900 Euclid Ave. 
Cleveland, OH 44106 

Eng Science Software Inc 
Attn: Dr. Kevin Walker 
P 65 Log Road 
Smithfield, RI 02917 

Exxon Res & Eng Co 
Attn: R P Gangloff 
Box 45 

Linden, NJ 07036 

FAA - New England Reg. 
Attn: Dan Salvano 
12 New England Exec. Park 
Burlington, MA 01803 

Failure Analysis Assoc. 
Attn: Dr. Jerrell Thomas 
2225 E. Bayshore Rd 
Palo Alto, CA 94303 

Grumman Aerospace 
Attn: R. Friedman 
Mail Stop B-43/35 
Bethpage, LI NY 11714 

General Electric Co 
Attn: Ted Russell G-25 
Turbine Division 
Schenectady, MY 12345 

NASA - LeRC 

Attn: J. Acurio, Director 

USAR&TL Propulsion Lab 
Cleveland, OH 44135 


Johns Hopkins Univ. 

Attn: Prof. W.N. Sharpe 
Chairman, Dept of Mech 
123 LaTrobe Hall 
Baltimore, MD 21218 

Lockheed Georgia 
Attn: H.S. Sweet 
Marietta, GA 30060 

McDonnell Douglas 
Attn: F.C. Claser 
Dept 243 
Box 516 

St Louis, MO 63166 

McDonnell Douglas 
Attn: R.T. Kawai (36-41) 
3852 Lakewood Blvd 
Long Beach, CA 90801 

Univ. of Minnesota 
Attn: Dr. W.W. Gerberlch 
Dept Chem Eng & Matl Scl 
Minneapolis, MN 

NASA - LeRC 
Attn: R.E. Gaugler 
Mail Stop 6-10 
Cleveland, OH 44135 

NASA - LeRC 
Attn: T.W. Orange 
Mail Stop 49-7 
Cleveland, OH 44135 

NASA - LeRC 
Attn: J.F. Saltsman 
Mail Stop 49-7 
Cleveland, OH 44135 

NASA - LeRC 
Attn: R.J. Simoneau 
Mail Stop 6-10 
Cleveland, OH 44135 

Northwestern Univ 
Attn: Prof S. Nemat-Nasser 
Dept of Civil Eng 
Evanston, IL 60201 

NASA - LeRC 

Attn: S&MT Div. Files 

MS 49-6 

Cleveland, OH 44135 


144 



DISTRIBUTION LIST (Continued) 


0 R M L Box X 

Attn: J. McGowan 
Oak Ridge, TN 37830 

Rensselaer Polytech Inst 
Attn: Dr. Erhard Krempl 
Mech & Aero Eng Dept 
Troy, NY 12181 

Rocketdyne 

Attn: Dr. D. Matejczk/ACIO 
6633 Canoga Ave 
Canoga Pk, CA 91304 

Rockwell Int'l 
Attn: Dr. Mike Mitchell 
Rockwell Science Ctr 
P.0. Box 1085 
Thousand Oaks, CA 91360 

Rocketdyne 

Attn: Jim Newell AC10 
6633 Canoga Ave 
Canoga Pk, CA 91304 

Sandia Labs 
Attn: D.W. Lobitz 
Div 1541 
Box 5800 

Albuquerque, NM 87115 

Southwest Research Inst 
Attn: Dr. Thomas Cruse 
PO Drawer 28510 
San Antonio, TX 78284 

Southwest Research 
Attn: Dr. C.S. Chan 
PO Drawer 285 
San Antonio, TX 78284 

Texas A & M University 
Attn: Dr. W.L. Bradley 
College Station, TX 77843 


Westinghouse - Hanford 
Attn: L.K. Severud 
Mgr, Plant Analysis 
Box 1970 

Richland, WA 99352 

Cummins Engine 
Attn: Justin Drake 
Tech Ctr 
500 S Poplar St 
Columbus, IN 47201 

University of Dayton 
Attn: Dr. Joseph Gallagher 
Research Institute 
Rm 563 Kettering Bldg 
Dayton, OH 45469 

Defense Documentation Ctr 
Cameron Station 
5010 Duke St 
Alexandria, VA 22314 

NASA-AFWAL Tech Liaison 
Attn: Everett Bailey 
AFWAL/DO 

Wriaht Patterson Air Force Base, uh 
a 45433 


Scientific Info Service 
Attn: Ref. USZ # 

7 Woodland Ave 
Larchmont, NY 10538 


0 R N L Box X 

Attn: Dr. Robert Swindeman 
Metals & Ceramics Div 
Oak Ridge, TN 37830 


NASA - LeRC 

Attn: N. T. Saunders 

Mail Stop 3-8 

Cleveland, OH 44135 


Texas A & M University 
Attn: Dr. Walter Haisler 
Aerospace Eng Dept 
College Station, TX 77843 

Williams Research 
Attn: R.L. Maxwell 
2280 W Maple Rd 
Walled Lake, MI 48088 


Scientific Info. Service 
Attn: Ref USZ 

7 Woodland Ave. 

Larchmont, NY 10538 


145 



GOVERNMt 


X RECIPIENT'S CAI 


arileWM 


4 TITLE AMD SUBlilU ” „ 

Creep Fatigue Life Prediction for Engine Hot Section 

Materials (Isotropic) 


" 7. AUTHORISE , . e 

Vito Moreno, David Nissley, Li-Sen (Jim) Lin 

9 PERFORMING ORG. NAML AND ADDRtbb 
UNITED TECHNOLOGIES CORPORATION 
Pratt & Whitney 
Engineering Division 

12. SPONSORING AGENCY NAME AND AUDKLbb 
National Aeronautics and Space Administration 
Lewis Research Center M1 « 

21000 Brookpark Road, Cleveland, Ohio 4413b 


i. REPORI uA it 
March 1 985 


PERFORMING ORG. REM. nu~ 
PWA-5894-34 

. WORK UNI I NO. 

. CONTRACT OR GRAN I NO. 
NAS3-23288 

TYPE REPr. /PERIOD CUVtKtb 
Annual Report (Second) 


P roject Manager^G^R . Halford, NASA Lewis Research Center, Cleveland, Ohio 


turbine hot section alloy (619UU Htj, lnaHino conditions and to define data 

crack initiation life for relevant creep-fatigue lQading conditions^an^^ ^^ ^ ^ 

required for oetermination of nodel const .ev the nodels satisfactorily 

rate, mean strain, strain hold times and teinperat . irements A cycle Damage 

:h t ;^ t I u ;"^d f Sllero:o e bse?»e;’ 0n 

SJKSS1. 

5,832 rrjr&ft 1-fl * 

The CUA ^el accurately predicts all of th t Trfdeteminable fro™ rapid 

cycled fuliy'reverse^fatigue 0 ^?' aSd nonotonic tensile and/or creep data. 


— 17 . KEY WORDS (SUGGESTED BY AUTHOR(s] 

creep fatigue , crack initiation 
life prediction. Isotropic material 


i. distribu i loH staiemtt 


19 . SCCURlIY CLASS THIS (tthPi mo. SCCUKlIV CLASS IHIG ll-AGEipi. HU. F 55 


:. PRICE *' 


* for sale by the National Technical Information Service, Springfield, VA 22161 


:|i ?t! I| Ill III 1 1 |i*l if 'f'l li|il| j, 4i!i 1 1 'I i* I'll" Mill lull !'i|illliil| Ijllii-l-lll-h-Hl i-IWltilHnH-i- in •( i i*:mn i-»h > I'i'iiln hi i '-i '■* i 'il l If III I i||i"|l- 1| l!j|i i |-|H if> « il-#* • HI jvlilN H-Ml-IrtiMi -I*! 1 • mu- « ■ illFkM't-i.i! -Mfi »iNmh»ii*'!MM-«H"I‘:*!||I-' *IN«lliWU|j-i-#(f#» i-i'IMk hi u 'Fllilhi'i' 





